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I. INTRODUCTION 


This report summarizes all work which has been performed by 
the Southwest Research Institute during the past two-year period under 
Contract NAS8-21133, n Study of Minimum Pressure Loss in High 
Velocity Duct Systems 11 . This study is being performed for the George 
C. Marshall Space Flight Center of the National Aeronautics and Space 
Administration and administered technically by the Propulsion and 
Vehicle Engineering Laboratory with Mr. R. H, Veitch serving as 
Technical Manager. 

The general objective of this study is to perform a theoretical 
and experimental investigation of ducting and ducting components to 
minimize pressure drop, surge pressures and vibration levels and, thus, 
to generally improve flow conditions of gases and liquids in space vehicle 
feed systems. Because of NASA needs, much emphasis has been placed 
on the study of bellows flow-induced vibrations. 

Summary of Result s 

A number of significant findings have been made throughout this 
study and are, discussed in detail in this report; these are summarized 
below: 

(a) Bellows P low-Induced Vibrations - The fluid -elastic 
mechanism causing bellows flow excitation (vortex shedding) 
has been observed and described. Analytical models have 
been developed to allow a designer to predict when flow 
excitation may occur, and to estimate the severity of the 
bellows vibrations. All pertinent data, in the form of 
equations and curves, presented in this report. The 
method is subject to certain current limitations which are 
discussed . 

(b) Bellows Liner Design - A limited amount of information 
is presented to aid in the design of a conventional bellows 
liner which suppresses flow-induced vibrations. 

(c) Bellows Pressure Loss - All available data has been 
compiled which gives bellows pressure loss for various 
convolution geometries, various sizes, and various flow 
media. Existing bellows pressure loss correlation methods 
have been reviewed, and while none are completely adequate, 
one is recommended for future use. 


(d) Elbow Pressure Loss - A new elbow design has been 
discovered which results in a significant reduction in 
pressure loss, as confirmed by tests. \ 

(e) External Damping Devices - Various bellows external 
damping devices have been tried as a means of suppressing 
bellows flow-induced vibrations. The results of these tests 
are presented to guide the designer in achieving vibration 
suppression by this means. 



II. BELLOWS FLOW-INDUCED VIBRATIONS 
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II. 1 Introduction 

The occurrence of flow-induced vibrations of metal bellows 
contained in fluid ducting systems has, for some time, been a problem 
for aerospace applications. There are known instances where flow- 
induced vibrations of bellows have resulted in fatigue failures which 
forced the premature shutdown of some critical fluid systems. The most- 
common methods for correcting unsuccessful bellows installations have 
been either to install an internal liner, where possible, or to use multiple 
plys or thicknesses of metal when constructing the bellows. Unfortunately, 
increasing the number of metal plys has not always cured the flow-induced 
vibration problem, and the use of a bellow's liner generally leads to an 
ii crease in component weight and cost. A major obstacle in the past has 
been that the flow mechanism which causes bellows vibration has not been 
described or understood so that flow-induced failures could be anticipated. 
From the standpoint of the designer, a desirable goal is to have available 
an analytical procedure which allows a prediction of critical fluid flow 
ranges for a givep bellows configuration and, further, which gives a method 
for estimating stress levels resulting from the flow excitation in these 
critical ranges. 

This chapter of the report is devoted to a comprehensive discussion 
of the bellows flow-induced vibration problem. As will be demonstrated, a 
fluid-elastic instability — -coupling of vortex shedding from the bellows con- 
volutions with the bellows structural vibration modes — is responsible for 
flow-induced failure problems. The conditions under which flow excitation 
is likely to occur are predictable, with acceptable accuracy, through the use 
of a simple analytical procedure. Further, good progress has been made in 
the development of a method for predicting magnitudes of the dynamic stresses 
which occur in the convolutions because of flew excitation. This refined stress 
prediction method is supplemented by a simple "stress indicator" procedure 
for quickly judging potential severity of bellows vibrations. A preliminary 
"stress indicator" versus number-of-cycles-to-failure curve rias been com- 
piled as a further aid to the engineer. Collectively, these various procedures 
represent a design or evaluation method which can be used to determine the 
adequacy of a given bellows application. 

II . 2 Background Information and Preliminary Test Results 

A search of the open literature, conducted at the start of the project, 
revealed that no previous work had been reported on flow excitation of bellows. 
In fact, the only discussions found in the literature which deal with the related 


problem of mechanically excited vibration characteristics of bellows are 
those of Daniels (1) and Lytle (2), A study of unpublished reports from 
the aerospace industry, however, revealed that bellows flow-induced 
fatigue failures had been observed for a variety of conditions, including 
liquid and gas flow, single - and multiple -ply bellows, and both free and 
wire braid -covered items. 

In an effort to better understand the problem, several different 
bellows specimens were initially tested in a water-flow loop. Figure 1 
shows a typical item mounted in the flow loop test section — in this case 
a free bellows with an internal diameter of 2-in. Each test bellows had 
one or more strain gages cemented on the bellows convolution crowns 
(see Figure 2 for bellows nomenclature) to allow flow-induced strain to 
be monitored. For each test bellows it was generally found that as the 
fluid velocity was slowly increased, starting from zero, successive longi- 
tudinal accordian vibration modes were excited. Figures 3 and 4 show' 
typical strain data and Figure 5 illustrates the vibration modes observed. 

As may be seen from Figures 3 and 4, each vibration mode was excited 
over some well defined fluid velocity range with the peak excitation occur- 
ring near the center of this range. A Fourier analysis of the strain signal 
was made at each test point with a wave analyzer, which revealed that, in 
some instances, more than one mode was being excited at a given fluid 
velocity. The simultaneous existence of two modes has only been observed, 
however, where the resonant frequency of one mode is about twice the 
resonant frequency of the other. In Figure 4, for example, the second 
mode showed significant response over the same fluid velocity rsmge for 
which the first mode was responding. 

While the brief description above is not complete as to many fine 
points of other observations about bellows flow-induced vibrations, it should, 
at least, introduce the reader to the problem. , 

II. 3 Flow Excitation Mechanism 


Based on test data of the type discussed above, and on the results 
of some flow visualization experiments, evidence is now conclusive that 
vortex shedding from the bellows convolutions is the prime mechanism 
responsible for the flow excitation phenomena which are observed in practice 
In every sense, the phenomena is a fluid -elastic instability (see References 
3-7), since a mutual interaction or coupling occurs between the dynamics 
of the bellows (the elastic system) and the nonsteady flow near and between 
the convolutions (the dynamic fluid system). Many factors have an influence 
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Figure 1. installation Of Vibration Test Beilows In Flow Loop 
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A mass < m ) is assigned each convolution crown 
and root; the number of masses is 2N C - 1 . 

The value of m is m m + mf where 

m m a 7rp m tNpD m 1 7T8 + ( h “ 2a ) I 


The number of springs is 2N C and 
k B 2N C k A 

where k A is the overall bellows spring rate. 


Figure 2. Convolution Nomenclature And Mechanical Vibration Model 
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Figure 3. Bellows Stress As A Function Of Flow Rate For Four 
Flow -Excited Modes Of Vibration - PN 08046 
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on this coupling phenomena; these factors will be discussed in detail in 
later sections of this report. The purpose of this section is to describe 
the fluid and structural behavior when vibration occurs, and also for 
those conditions when it does not occur. 

For flow through a convoluted section of bellows, with the convo- 
lutions fixed to prevent vibration, a steady turbulent flow situation exists 
throughout the bellows. Figure 6 illustrates the steady flow in a bellows 
with no convolution vibration; this is a physical picture formulated from 
our own flow visualization experiments, and from the information reported 
by Haugen and Dhanak {8), and Townes and Sabersky (9). Each convolution, 
in general, experiences an internal eddying flow because of a momentum 
transfer process. Energy is transferred into each convolution in a mixing 
zone or shear layer, indicated in Figure 6. The number and size of the steady 
vortices which exist in the convolutions are a function of the convolution 
geometry, also indicated in Figure 6. Very shallow convolutions may con- 
tain only one vortex, while two or more may exist in deep convolutions. As 
will be discussed in a later section of this report, it is this momentum 
transfer process and vortex motion which causes a bellows to have an extra- 
ordinarily high pressure loss. 

The physical picture described above is strictly for the case of no 
convolution motion or vibration. The only flow excitation which the bellows 
would experience with this type of flow behavior, would be random as a 
result of the turbulent pressure fluctuations experienced at the fluid - 
convolution boundary. This excitation mechanism would give forces of a 
much lower magnitude than the one really responsible for bellows flow 
excitation. 

When the bellows is allowed to be flow excited, the physical picture 
of internal flow described above is no longer valid, in general. For this 
case, the flow now involves a large-scale periodic vortex formation and 
shedding process which has been visualized With the apparatus shown in 
Figure 7. The model consists of a two-dimensionai clear plastic channel 
containing a .short section of convoluted metal which is placed so that the 
convolution tips are exposed to the fluid flow, as in an actual bellows. One 
end of the convoluted segment is movable, in a longitudinal direction, so 
that the convolution pitch can be'varied, as desired, by stretching or com- 
pressing the segment. When water is passed through the channel at the 
proper velocity, the segment is flow excited (if the conditions are right); 
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Figure 6. Fluid Flow Behavior With Stationary Convolutions 
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this vibration of the convoluted segment is easily viewed in a slowed - 
down manner with the aid of a strobe light. By injecting ink upstream 
from the convoluted segment, as indicated in Figure 7, the vortex 
shedding process can readily be visualized. 

Both still and moving pictures of the vortex shedding have been 
made; Figure 8 shows several of the still pictures, which were taken at 
random times during the vortex shedding process. The presence of a 
vortex pattern is evident ill each picture; The movie, of course, shows 
the interaction of the fluid motion and the segment much more clearly 
than do the still pictures. 

Figure 9 shows the sequence of fluid and convolution events which 
have been observed in a frame-by-frame examination of the motion 
picture. Note, first of all, that the mode of vibration of the segment is 
one where each convolution moves out-of-phase with the adjacent convo- 
lutions). The vortex shedding process on the vibrating segment, as 
shown by Figure 9, occurs as follows: 

Postion I - A large vortex (c) has formed between convolutions 
1 and 2, and is being "pushed" out into the fluid stream by the 
pinching action of these two convolutions. A large, well-formed 
vortex (a) is moving across the tip convolution 3; the origin of 
this vortex will be made clear in the description of other events. 
Finally, a small vortex (b) is beginning to form on the downstream 
side of convolution 2. 

Postion II - Vortex (c) has been pinched from between convo- 
lutions 1 and 2 and has now moved out into the fluid stream, and 
will soon "detach" from convolution 1 and be swept downstream. 
Vortex (b), which is forming on the downstream side of convolution 
2, is gaining in strength. Vortex (a), which in Position I was 
beginning to move over the tip of convolution 3, is now further 
downstream. 

Position III - The space between convolutions 1 and 2 is opening 
tip, hence fluid is being drawn in, and vortex (d) is beginning to 
form on the downstream side of convolution 1. Vortex (c) has been 
swept over the top of convolution 2, with the timing being about 
right to reinforce vortex (b), which is being pinched from between 
convolutions 2 and 3; the combined vortex is labeled (b + c). 
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Figure 8. Series Of Photographs Showing Vortex Shedding From 
Two-Dimensional Convoluted Metal Segment 





jj Figure 9. Sequence Of Coupled Fluid*- Convolution Events Observed 

With Two-Dimensional Beilows Flow Visualization Model 
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Position IV - Vortex (d) has gained in strength as it forms on 
the downstream side of convolution 1. Vortex (b + c) is being 
pinched out from between convolutions 2 and 3, and will begin 
to propagate downstream. 

Position V - This position is the same as Position I, so the 
cycle is complete. Note that vortex (a) in Position I was a 
combination of two vortices, as is vortex (b + c) in Position V. 

It must be emphasized that the series of fluid events described in 
the foregoing are probably not exactly descriptive of other situations. 
Probably the main difference seen from a comparison with other situations, 
would be that the phasing of vortices propagating downstream would not 
always be correct, with respect to a vortex forming at a point downstream, 
to give a reinforcement condition, as was seen here. We have some evi- 
dence to suggest that there might be optimum convolution geometries, for 
a given mode of vibration, which yields a peak excitation amplitude. This 
might be caused by a vortex reinforcement phenomena. 

It has been shown in the foregoing discussion that two different 
fluid flow regime'? exist; one occurs when the convolutions are constrained 
and, consequently, no organized large scale vortex shedding is seen, only 
a steady-state turbulent situation. The other flow regime occurs when the 
convolutions are vibrating as a result of flow excitation; here a large- 
scale vortex shedding process does exist. Therefore, the bellows flow 
excitation can be described as a mutual or coupled instability between the 
bellows structure and the fluid vortex shedding process. One cannot exist 
without the other. 

The primary condition which must be satisfied for this coupling to 
occur is that the bellows structural frequency and the ideal vortex shedding 
frequency be equal, o:d about equal within certain limits. This limit of 
vortex-structure coupling can be called a "lock-in" range, and as shown in 
Figures 3 and 4, it can exist over a rather broad velocity range. 

II. 4 Prediction of Flow-Excitation Lock-In Range 

The purpose of this section is to discuss how the occurrence of 
the vortex shedding and bellows structure coupling can be predicted. 

Basically, this involves calculating the ideal vortex shedding frequency 
and he bellows vibration mode frequencies, and then comparing the two. 

Also, some allowance must be made for the fact that a "lock-in" range exists. 
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Vortex Shedding Frequency 


For a given bellows, having some particular convolution geometry 
and structural mode frequencies, there are certain optimum fluid velocities 
which result in a maximum amplitude bellows excitation of each mode (see 
again Figure 3). It is at these velocities that the vortex shedding process 
is best able to feed energy from the fluid stream into the vibration process. 
In other words, the vibration frequency and fluid velocity conditions are 
"optimum" from a vortex shedding standpoint. It has been found that the 
use of a Strouhal number is an excellent means of correlating the vibration 
frequency, fluid velocity and geometry under these optimum conditions, as 
is true for any vibration phenomena involving vortex shedding. 

In general terms, a Strouhal number is a dimensionless quality of 
the form 


S = f 4/V (1) 

where f is a frequency, I is a length quantity, and V a fluid velocity. 

For the case of bellows flow-induced vibrations, the only problem in using \ 
this correlation parameter is in selecting a satisfactory length quantity. 
Referring back to the discussion of the vortex formation and shedding process 
illustrated in Figures 8 and 9, it would seem that two length quantities may 
be primarily involved; these are the convolution pitch and the convolution tip 
width a . While neither of these quantities has been found entirely 
satisfactory, that is neither gives a precisely constant value of S for all 
conditions, we have chosen to use a in making Strouhal number calculations. 

Figures 10 and 11 show example data of the Strouhal number calcu- 
lated at the optimum flow excitation condition for two different bellows. 

The pitch of these bellows were changed by stretching and compressing them 
from their original configuration. Strouhal number values based on both 
the tip width a and the pitch \ have been plotted. Note how neither of 
the two dimensionless numbers is constant as pitch changes. Based on a 
composite of all test data, a single curve of Strouhal number (based on tip 
width a ) as a function of the ratio of pitch to tip width ( \/o ) has been 
prepared and is believed valid for a general case. This curve is shown in 
Figure 12; also shown are dotted curves defining the limits of the lock-in 
range. Using this data, it is possible to calculate the fluid ve 1 deity V at 
which a given bellows longitudinal vibration mode (defined by a frequency 
value) will most likely be excited and, also, to estimate the total lock-in 
range of a givv..n mode. To make this calculation it is necessary to know 
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Figure 11. Strouhal Number Correlation For Bellows #105 ( first mode ) 








three quantities; namely , the pitch \ , the convolution tip width a and 
the mode frequency f m . The value of this critical velocity is, then 

V = (2) 

s a 

where S a is obtained from Figure 12. 

Bellows Structural Frequencies 

In order to calculate the fluid velocity range(s) over which bellows 
excitation will occur by the procedure described in the foregoing section, 
it is necessary to know the frequencies of all bellows modes which may be 
flow excited. Therefore, some method must be used to calculate these 
frequencies. A method, ba.sed on an analogous mechanical model, has 
been developed and will be discussed below. Because the structural 
bellows modes cannot be considered independent of fluid loading, this 
influence is discussed also. 

It has been determined from experiments that three different kinds 
of structural modes can be flow excited, these are; axisymmetric longi- 
tudinal or accordiar modes, longitudinal antisymmetric or cocking modes, 
and local convolution bending modes. These different kinds of bellows 
modes are illustrated in Figure 13. The longitudinal modes (both sym- 
metric and antisymmetric) involve vibrations where there is longitudinal 
motion of the roots and/or crowns of the convolutions. The higher order 
local bending modes will involve no longitudinal motion of roots and crowns, 
only local flexural distortion of the convolutions. In our experiments, we 
have observed only one type of higher order bending mode; it is the one 
shown in Figure 13 and again in -Figure 14. It is doubtful that many, if any, 
other higher order bending modes can couple with vortex shedding to pro- 
duce a bellows flow-induced vibration; this has not been confirmed experi- 
mentally. For the purposes of this report, we will assume that the higher 
ordsr bending mode illustrated in Figure 14 is the highest which can be 
excited by vortex shedding. The reader should keep in mind, however, that 
this does not rule out the possibility of other higher order local convolution 
bending modes being excited. 

During the course of the present study, the mechanical model 
shown in Figure 2 has been found to yield adequate predictions of the longi- 
tudinal symmetric mode frequencies of a free bellows, if the overall spring 
rate is accurately known. The [model consists of (2N C ~1) elemental masses 
m and 2N C elemental springs k connected alternately in series, giving 
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Figure 14. Only Observed Higher Order Mode 
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an N = 2N C -1 degree-of-freedom system; N c is the number of bellows 
convolutions. Each elemental mass represents one-half of the mass of 
metal contained in each bellows convolution plus some fluid added mass, 
and k represents the spring rate of one -half of a convolution. A 
tabulation of dimensionless frequencies obtained for this model up to 
25 degrees of freedom is shown in Table I; the normalizing frequency is 

JL /JA 1/2 

° 2 it V m / 

For use in the mechanical model, the values of the elemental 
spring rate k and the elemental metal mass m m may be estimated by 
the following procedure. The value of k is defined in terms of the 
overall bellows spring rate Ka 

k = 2 N c K a (3) 

where K A is known either from a spring rate test or from some analytical 
calculation.^ The value of the elemental metal mass m m is equal to 
one-half of the mass of one complete convolution and may be satisfactorily 
approximated by the expression 


m m = "Pm 1 N p D m [' rra+(h " 2a !] ( 4 ) 

where (see Figure 2) 

mass density of bellows metal 
thickness per ply of metal 
number of bellows plys 

mean diameter of bellows convolution. D m =(D^+D 0 )/2. 
mean radius of convolution 
convolute height 

The frequencies obtained using the information given above are 
for the symmetric longitudinal modes only. An analysis has been made to 
determine, the cocking mode frequencies and the results show that these 



Appendix B gives a review of existing bellows spring rate information. 
In general, it has been found that experimentally obtained spring rate 
data is the most satisfactory. 
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new modes are always within the frequency limits of the symmetric 
longitudinal modes. Therefore, a calculation of the range of symmetric 
longitudinal modes can be considered to include any cocking modes which 
may occur. Cocking modes have been observed only for bellows with 
quite deep convolutions, and will probably not occur in most instances. 

The higher order local bending mode shown in Figure 14 may be 
modeled with the mechanical analog shown in the same figure. 

Fluid Loading Effects 

The previous discussion concerning calculation of structural 
bellows mode frequencies made no mention of possible fluid loading 
effects. The presence of the fluid can, however, cause some changes 
in the structural frequencies and even add more degrees of freedom to 
the system (increase the number of possible modes). 

When a bellows is being excited by a flowing fluid, two basic 
classes of forces are present; namely, active forces which are vortex 
shedding derived, and passive fluid forces which are present regardless 
of the excitation source. For the idealized case of an incompressible 
fluid, only two types of passive fluid forces are possible and these give 
rise to an added-mass and damping influence; t!\ey result because the 
fluid must be moved in and out from between the convolutions and, to 
some extent, back and forth longitudinally as the bellows vibrates. If 
only added-mass and damping effects are introduced, the number of 
degrees of freedom of the bellows (that is, the number of possible vibra- 
tion modes) does not necessarily change. One exception to this are in-phase 
and out-of-phase flex hose modes to be discussed later. 

For a compressible fluid, the complexity of the situation is greatly 
increased, because the bellows' vibrations may couple with various system 
acoustic modes u Generally, this increases the number of degrees of 
freedom, or the number of possible vibration modes of the bellows. This 
phenomena is illustrated in Figure 15 which shows the response character- 
istic of a single convolution test model as influenced by an acoustic loading 
effect. One curve in this i'igure shows the frequency response amplitude 
of the test convolute where no significant ducting longitudinal acoustic 
resonance occurred; the response is clean, and has a single peak at the , 
expected frequency of 600 cps. The mechanical analog for the convolute 
with no acoustic effect is also shown in Figure 15; additional mass and 
damping have been added to simulate these fluid loading effects, The other 
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Figure 15. Vibration Response Of Convolution Test Model Ring Number 2 In 
Water ( No Flow ). With And Without Ducting Acoustic Resonance - 
Ring Excited Electromag netica I ly 


27 




1 


curve in Figure 15 shows the response of the test ring where the effect 
of a duct acoustic resonance was significant. Note that two distinct 
resonant frequencies are present. In this case, a description of the 
convolute dynamics must include the line acoustic loading effect. While 
the system was ordinarily considered representable by a single-degree- 
of-freedom second order system, it must now be considered as at least 
a two-degree-of-freedom system; a possible mechanical analog is shown 
in Figure 15. 

For the present, only incompressible fluid loading effects will 
be considered in regard to coupled bellows-fluid mode frequency calcu- 
lations. Further, because the fluid damping is small enough to give a 
very small change in the system resonant frequencies, it will be dis- 
regarded in all frequency calculations. Later, however, damping will 
be shewn to be of great importance in calculating bellows flow-induced 
stress levels. 

For the first few modes of a free bellows (more precisely, where 
the mode number N is small with respect to the number of convolutions 
N c ), it has been determined by comparing calculated with observed 
frequencies that the fluid added mass is adequately represented by setting 
it equal to the total amount of fluid trapped between adja-cent convolutions; 
see Figure 2. Physically this makes sense since the first few modes 
involve a gross back-and -forth motion of adjacent convolutions with the 
trapped fluid carried along with the convolutions. This means, then, that 
each elemental mass of the mechanical model must include one-half of 
the mass of fluid trapped between two adjacent convolutions, or approximately 

m f = 2 P f D m h ( 2a " tN v ) (5) 

with, now, the elemental model mass becoming 

' \ 

m = m£ + m m ! (6) 

where the value of m m is the metal mass discussed previously. In 
Equation (5), Pf is the fluid mass density; other symbols were defined 
previously. 

For the higher modes of a free bellows, it is not adequate to 
consider the fluid added mass as being simply that amount of fluid trapped 
between convolutions. The fluid behavior for the higher modes is rather 
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more of a "squeezing in-and-out" from between the convolutions type of 
motion. This acceleration of the fluid out of and back into the convolutions 
results in alternating pressure forces which may be interpreted as an 
added-mass effect with respect to the bellows vibrations. In Appendix C 
this added mass is analytically estimated to be 

wD m<>f h3 

m f = ~6 < 7 > 


where 6 is the gap width between adjacent convolutions as shown in 
Figure 2. 

Comparison of Calculated and Observed Frequencies 

Comparison of observed and calculated frequencies for several 
different free bellows specimens has verified the validity of the free bellows 
frequency calculation model described in the foregoing. If the spring rate 
of a given bellows is accurately known from a force-deflection test, the 
calculated frequencies, where the fluid added mass is small, are surpris- 
ingly accurate; generally the error is no greater than 10 percent and often 
no more than 5 percent compared with experimental values. When the 
fluid added mass is of the same order of magnitude as the metal mass or 
greater, the errors begin to increase in some instances, particularly where 
the bellows geometry is unusual. For example, if the convolution is an 
extreme "open" design, the added fluid mass is overestimated by the 
./ ' equations given previously. This is understandable since the fluid cannot 
be readily "held" between the convolutions as they vibrate back and forth. 

Table II gives examples of experimental and calculated frequencies 
for many different cases. In most instances the accuracy obtained from 
using the mechanical model described above is completely adequate. More 
sophisticated computational methods could well be employed, however, 
because of uncertainties in various bellows dimensional data when comparing 
"blueprint" information and real bellows dimensions, very little, if any, 
additional accuracy can be achieved. The largest source of error in using 
the mechanical model for calculating frequencies, is an inaccurate knowledge 
of the bellows spring rate; see Appendix B v 

Flex Hose Modes 

The difference between the convolution vibrations of a flexible 
hose and a free bellows is that the allowed modes are limited by the wire 
braid covering. Inmost cases this covering clamps tightly around the 
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t 


Bellows 

Internal 


n 

Frequency (cps) 


Number 

Media 

Excitation | 

Calculated 

tsssssBSsi 

Error 








#101 

Water 


Flow 

504 

494 

2. 0% 

#104 

Water 


Flow 

770 

813 

5. 3% 

#105 

Water 


Flow 

527 

515 

2. 3% 

#106 

Air @ 30 psig 


Shaker 

790 

935 

15. 5% 

#106 

Water @ 30 psig 


Shaker 

616 

700 

12. 0% 

#106 

Water 


Flow 

616 

764 

19.4% 

#110 

Water 


Flow 

78o 

836 

6. 7% 

#112 

Air @ 30 psig 


Shaker 

635 

620 

2.4% 

#112 

Water @ 30 psig 


Shaker 

521 

461 

13. 0% 

#112 

Methanol @ 20 psig 

Shaker 

541 

478 

13. 0% 

#112 

Methanol @ 30 psig 

Shaker 

1065 

{2nd Mode) 

914 

16. 5% 

#112 

Water 


Flow 

521 

521 

— 

#114 

Water 


Flow 

.774 

754 

2.7% 

#114 

Water 


Flow 

1512 

1435 

5.4% 


Note: See Appendix A for Bellows Dimensions 


TABLE II '• '' 

Calculated and Experimental Frequencies 
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bellows when the hose is pressurized, preventing longitudinal motion of 
the convolution crowns. Except for special instances where not all con- 
volutions crowns are in contact with the wire braid, it is believed that the 
only purely structural longitudinal mode which can exist is the one shown 
in Figure 16(a). Because of fluid addect-mass effects, however, two 
coupled fluid-structural longitudinal modes appear possible, and these 
are shown in Figure 16(b) and (c). Note that botl:| of the coupled fluid 
structural modes are possible for the single allowed structural mode. 

The primary difference between the two coupled fluid -structural 
modes shown in (b) and (c) of Figure 16 is the fluid mass loading involved. 

For the case in Figure 16(b), the in-phase mode, the fluid-added mass 
corresponds to a portion of fluid contained between adjacent convolutions 
and may be estimated from Equation (4). This quantity of mass accounts 
for the fact that fluid is accelerated back and forth with the convolutions 
as they vibrate in the in-phase mode. For the out-of-phase mode, the fluid 
added mass does not correspond to some real quantity of fluid, but is an 
appareTit mass which resu?A horn the fluid being accelerated in and out 
from between the convolu * ' ■ \S as discussed in the foregoing for the free 
bellows. This quantity of added mass may be estimated from Equation (7). 

In addition to the longitudinal flex hose modes discussed above, it 
is believed possible for the higher order local bending mode shown in Figure 
14 to also exist. 

Thus far, only one flexible hose has undergone extensive flow testing. 
This particular item, a 1.5-in-id specimen, was placed in the water flow 
loop and subjected to a slowly increasing internal fluid velocity while bellows 
response was noted. Because the convolutions were covered on the outside 
with the wire braid, it was impractical to measure convolution strain, or 
even observe the convolution motion stroboscopically as had been done with 
the free bellows. The bellows response was therefore detected indirectly 
by three methods; external acoustic emission, dynamic fluid pressure down- 
stream, and acceleration of the duct adjacent to the hose were observed. The 
external acoustic emission technique gave the most reasonable results for 
this test but, of course, gave only a qualitative indication of response. 

Prior to the te^t, calculation of the in-phase and out-of- phase 
mode frequencies gave values of 6720 cps and 3420 cps, respectively. 

The measured spring rate was used in this frequency calculation. Notice 
that, by coincidence, the ratio of these two frequencies is about 2 to 1. 

All of the output data from each monitor was subjected to a Fourier analysis 
to allow the bellows response modes to be clearly seen. In general, two 
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Provided By Braid 
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( a ) Structural Mode For Flexible Hose 



Figure 16. Summary Of Flexible Hose Longitudinal Vibration Modes 
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flow-induced modes were observed, one at about 3300 cps and the 
other at twice this value or 6600 cps. Figure 17 shows a plot of 
typical spectral data obtained from the external acoustic monitor to 
show how these modes responded as the flow velocity was increased. 
Figure 18 shows a corresponding plot of acoustic emission level as a 
function of Strouhal number. 

The fact that two modes were observed in this test of a flexible 
hose, very close to the calculated values for the in-phase and out-of- 
phase modes, substantiates to some extent the existence of these two 
types of modes. The fact that one happened to have twice the frequency 
of the other was pure coincidence; they were both apparently excited 
simultaneously because of the multiple -frequency content of the vertex 
shedding. 


Estimating Flow Excitation Range 

The foregoing discussions have described how the conditions for 
bellows excitation can be correlated with a Strouhal number, and, further, 
how the modal frequencies of a bellows can be calculated. The purpose of 
this section is to summarize how this information can be used to predict 
possible bellows flow-excitation ranges. 


As we have seen, a given bellows will have a number of structural 
modes which may be flow excited. These modes may be defined in terms 
of a number of modal frequencies denoted by 

1 * ^m3 * • - * fmn 


Here f m i denotes the first or lowest modal frequency and f mn denotes the 
highest. For the case of a free bellows, f m ^ will represent the first 
longitudinal mode and f mn will denote the convolute bulging mode; the 
total number of modes will depend on the number of convolutions. For a 
flex hose, there are, so far as we presently, know, only three possible 
modes; these are the in-phase and out-of-phase longitudinal modes shown 
in Figure 16, and the convolution bulging mode shown in Figure 14. Table 
III and IV summarize the frequency calculation procedure for free bellows 
and flex hose. 


Each of the bellows modes may experience flow excitation over a 
fluid velocity range defined by 

f m v 


V, = 


V = 
u S 


^cru 


oi 


( 8 ) 

( 9 ) 
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Acoustic Level -Arbitrary Units 







Step A 


Consider the bellows structure representable by the 
lumped mass-spring mechanical model shown in Figure 2, 


Calculate the elemental spring rate value k from the 
expression 

k = 2 N c K a 

where is the overall spring rate determined from 

a method given in Appendix A or from a force-deflection test. 

Calculate the elemental metal mass m m from Equation (4) or 

m m = 1T P m t N P D m [" a + ( h -2a)] 

Calculate the fluid added mass mf, using Equation (5) 

• for the first few longitudinal modes (N values) and 
Equation (7) for the higher longitudinal modes; thus 

m f = ~ p£ D m h (2a-tNp)J First few N values 


ir D m p f h : 


Higher N values 


Calculate the reference frequency f Q from the expression 

f = jl jjt 

° 2 tt V m 


where 


rn = m m + mf 


Look up proper dimensionless frequencies in Table I 
(corresponding to given N c value) and then calculate the 
tree mode frequencies by multiplying by the reference 
frequency. This gives all of the longitudinal mode frequencies, 


Calculate the first bending mode frequency from the informa- 
tion in Figure 14, or 

f s _L Ip? = _L JW 

2tt V mT" 2tt V m 


where 


k = 2k 


and m m +m£ 

m' - m/2 = "~ m ~2 

For the value of m£ use 

mf = irD^pfh^/36 

i ^ mh > ■iii h 'iii h . i b ii ii f i i i M i i ii ui.'n ' i nn " ii ii i ■ 11 — “ — ^ ~ " 

TABLE III - SUMMARY OF FREE BELLOWS FREQUENCY 
CALCULATION PROCEDURE (ASSUMED 
INCOMPRESSIBLE FLUID) 
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Step A 


Step B 


Step C 


Step D 


Step E 


Step F 


Consider the hose convoluted structure representable 
by the mass-spring model shown in Figure 16. 

Calculate the elemental spring rate value k from the 
expression 


k x 2 K a 

where is the overall spring rate of one complete 
convolution. The value of may be determined by a 
method given in Appendix A or from a force-deflection test. 

Calculate the elemental metal mass m from the expression 

1 m m — '"Pm*' 1’ (h— 2a)] 

Calculate the in-phase mode and out-of-phase longitudinal 
mode fluid masses from 

m f = \ PfD m k (2a-tNp)j in-phase mode 

vDmPfh 3 \ 

m f ” j out-of-phase mode. 

Calculate the in-phase and out-of-phase longitudinal 
frequencies from 



v/here 


m = m m + m£ 

Calculate the first bending mode frequency from the 
information in Figu re 1 4, or 

JL F _ j_ g 

2ir V m' ™ 2tt V m 
k' = 2k 

( m m m ^ m f 

m = 2 = 2 

For the value of mf use 

m f ” frD m p f h 3 /36 


where 


and 


TABLE IV - SUMMARY OF FLEX HOSE FREQUENCY 

CALCULATION PROCEDURE 
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and the optimum or most severe excitation will occur at a velocity 
equal to 


fm * 7 


V = 


In Equations (8) to (10), we have 


m 


S 


cru 


oi 


a 


V, 


V. 


u 


( 10 ) 


a bellows mode frequency 
upper limit Strouhal number 
lower limit Strouhal number 
optimum Strouhal number J 
velocity at lower limit of lock-in tfange 
velocity at upper limit of lock-in range 
convolution tip width 


see Figure 12 


Because of the lock-in range phenomena, and the closeness of successive 
b ellows mode s , it is best to assume that flow excitation can occur con - 
tinuously over a veloc it y range from the minimum velocity for the lowest 
mode to the maximum velocity for the highest mode . Therefore, if f m j 
and f mn are the lowest and the highest mode frequencies, the extreme 
possible limit of flow excitation is defined by the limiting velocities 


and 


V. = 


f ml c 


cru 


( 11 ) 


V 




u 


S 


a.i 


( 12 ) 


In summary, the range of possible bellows flow excitation may be 
predicted as follows: 


(a) Calculate the lowest and highest bellows mode frequencies, 
f , and f . The procedure for doing this is summarized 
in Tables in and IV. 


(b) Calculate the limits of fluid velocity corresponding to these 
two frequencies. This is done with Equations (11) and (J2), 
and with the values of the Strouhal numbers S aJ g and S aU 
obtained from Figure 12. 
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(c) 


Compare this flow excitation fluid velocity range with 
the known operating range of the bellows. If an overlap 
of these ranges exists, then excitation may occur. 


A graphical method of predicting bellows excitations ranges is 
illustrated in Figure 19. This method involves preparing a plot of fre- 
quency versus fluid velocity which contains both the bellows mode 
frequency information and the possible vortex shedding frequency limits. 
The upper and lower bellows frequencies are represented by horizontal 
lines on the plot (constant frequency lines), as shown in Figure 19. Also, 
lines of constant slope which pass through the origin are drawn :.n to 
represent the vortex frequency limits (as a function of velocity), defined 
by the expressions 



S of V 


<y 


(13) 


and 


f u = 


S V 

cru 


(14) 


The range of possible flow excitation is clearly seen from this plot. 
II. 5 General Discussion of Bellows Forced Vibration 


All previous discussion has been directed to the problem of determine 
ing when a bellows may be flow excited, based on a coincidence of coupled 
bellows -liquid frequencies and the vortex shedding frequency. This coinci- 
dence of frequencies does not guarantee that the bellows will vibrate, however, 
or indicate the' severity of vibration. It is desirable, therefore, to be able 
to predict vibration amplitudes and corresponding stress levels for a given 
bellows and given flow conditions. 

Figure 20 illustrates the fluid and structural dynamics involved in 
bellows flow excitation. The process of periodic vortex formation and shed- 
ding causes a corresponding periodic pressure to be exerted on the convolu- 
tions. The amplitude of this alternating pressure is proportional to the 
free-stream stagnation pressure (1/2 p^V^), So far as the bellows structure 

is concerned, the effect of this alternating pressure may be considered as a 
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Frequency 



Fluid Velocity 

2398 


Figure 19. Frequency versus Velocity Plot To Illustrate 
Bellows Flow Excitation Region 



net force applied at the tip of each convolution, as shown in Figure 20. 
The amplitude of this force is assumed to be of the form 


F = CpA p (1/2 p f V 2 ) 


(15) 


where Cp is a vortex force coefficient (a dimensionless coefficient) and 
A p is the projected convolute height area* over which the pressure acts 
(the fluctuating pressure producing the force). 


If the frequency of the vortex shedding coincides with some bellows 
mode frequency, then a resonance will occur and the bellows convolutions 
will experience a vibratory displacement, as illustrated in Figure 20, of 
the form 


x 


Cm F Q 


k a 


( 16 ) 


In Equation (16), C m is a dimensionless factor dependent on the~bellows 
mode of vibration (mode factor), F is the vortex shedding force defined 
in Equation (15), is the bellows overall spring rate, and Q is the 

dynamic amplification factor (damping). This convolution vibratory dis- 
placement x will cause a corresponding stress, and the magnitude of 
this stress is dependent on the convolution geometry, as illustrated in 
Figure 20. 

\ 

Accepting this simple concept of bellows vibration resulting from a 
coupling with the vortex shedding, the problem of calculating resultant 
stress levels boils down to one of requiring a forehand knowledge of the 
various "factors’ 1 illustrated in Figure 20. The primary factors are: Cp, 
a vortex force coefficient; C m a vibration mode factor; Q a dynamic 
amplification factor; and C s , a geometric stress factor**. Other factors 
may also be introduced to account for various unique situations. For 
example, when an elbow is located upstream of a bellows, this can have 
an effect and cause higher- than -normal stress levels for a given flow 


In earlier reports we used the convolution projected tip area as the 
area over which the vortex shedding fluctuating pressures act; we have 
since found that the convoluted projected height area is more realistic. 

## 

The stress model indicated in Figure 20 results from the analytical 
form presented by Salzmann (10), and is only representative of a 
number of stress prediction methods (see Appendix B). 



Vortex Shedding From Convolutions 



In the above equations 

Cp * vortex force coefficient 
Cp a elbow factor 
C m » vibration mode factor 
Q = dynamic amplification ( damping ) 
C s - geometric stress factor 


Figure 20. Illustration Of Stress Resulting From Vortex Force 



velocity. This can be accounted for by an "elbow factor" Cg, as 
indicated in Figure 20. If an acoustic resonance exists, this can change 
the picture, also. The objective of the next several sections is to give 
detailed results of studies of these various factors. 

II. 6 Vortex Force Coefficient - Idealized Test Model 

The initial study of the nature and magnitude of the vortex shedding 
force was undertaken with a special single convolution test model which 
is illustrated in Figure 21. A bellows convolution is simulated in this test 
model by a flexible ring which is geometrically quite similar to a convo- 
lution. The ring was constructed from steel and had a coating of lead on 
the flexible web section to provide sufficient damping to prevent early 
fatigue failure during actual testing; this happened with the first several 
rings w-hich were constructed without added damping. The test ring or 
simulated convolution is clamped in a special housing between a pair of 
exciter coils; also, a displacement probe is built into the apparatus to 
allow ring vibration amplitudes to be monitored. 

With this apparatus the following types of tests have been performed: 

(a) frequency response of the ring in air with excitation provided by the 
coils; (b) freouency response of the ring in water with excitation provided 
by the coils; and (c) response of the ring as a result of flow excitation wis.h 
internal water flow. Also, a calibration of static force versus deflection 
was made for the ring. From these tests and the force -deflection calibra- 
tion, several kinds of information were obtained. First, the forced vibra- 
tion test in air gave essentially the ring-only natural frequency and damping 
(neglecting air loading). Second, the forced vibration test in water allowed 
the water added-mass and water damping to be calculated. Finally, the 
flow-induced vibration data allowed an effective vortex force coefficient to 
be obtained over the lock-in range. This effective vortex coefficient was 
reduced from the ring flow-induced vibration data in the following manner: 

(a) The vibration amplitude (as a function of fluid velocity) was 
converted to an apparent force through the use of the force- 
deflection calibration. 

(b) This apparent force was converted to a true vortex force 

>y dividing by the Q value reduced from the forced vibration 
test data (in water). 

(c) This true force was converted to a force coefficient through 
the use of Equation (15). 
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x(t) 

SINGLE BELLOWS CONVOLUTION 



MECHANICAL MODEL OF CONVOLUTION 



Figure 21. Diagram Of Convolution Vibration Test Model And Equivalent 

Bellow^ And Mechanical Model Representations 
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Vortex coefficient data was obtained in this manner for two different 
convolution geometries and the results are shown in Figure 22. Notice 
the great reduction in the force coefficient as the convolution pitch is 
"opened up". 

II. 7 Vibrations of Real Bellows 


Bellows Mode Factor 

Having available the results of the experiments with the single 
convolution test model, discussed in the previous section, the next logical 
step was to use this force coefficient information, along with the equiva- 
lent mechanical model concept, to predict flow-induced vibration amplitudes 
of real bellows. This was done by assigning one force per mass point*, as 
shown in Figure 23, and then calculating the forced amplitude of the mech- 
anical model. Comparing the resultant vibration amplitudes with observed 
values showed the predictions to be high by a factor about equal to the 
number of convolutions involved in each "half wavelength" of bellows 
vibration. This suggested, therefore, that there was only one effective 
vortex force per mode number of the vibration. Further detailed com- 
parison of calculated and observed vibrations of bellows confirmed this 
idea. Also, reexamination of the visualized vortex shedding activity 
(illustrated typically in Figure 8) and reconsideration of the pressure 
forces involved, showed how this must happen. 

Figure 24 depicts the action of a bellows with several convolutions 
vibrating in the first longitudinal mode. As shown, the vortex shedding 
causes fluctuating pressure forces on each side of each convolution. 

Because of the phasing of these pressure forces, '^owever, only one 
effective vortex force exists per mode number of vibration. The first 
mode has one effective force point, the second mode, two, et- . Each 
of these effective forces is exerted a* the central convolution of each mode 
half wavelength, or at the point(s) of maximum displacement of each mode. 


By assigning one force per mass point, where each mass point 
represented one-half of a convolution, it was necessary to take 
each force equal to one -half of the force in Equation (15). 


/ 
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Vortex Force Coefficient - 


i 

t . 



Strouhal Number - fa/ V ™ 

Figure 22. Vortex Shedding Force Coefficients From Tests Of 
Two Simulated Convolution Geometries 
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All forces are of the form f n a f o sin(«t + 0 n ) 
where p depends on mode of vibration 



j Figure 23. Mechanical Model With One Force Per Mass Point - 

* • All Forces Are Of Equal Magnitude 

) 
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Equivalent Mechanical Model 


2353 


Figure 24. Effective Vortex Force On Bellows Vibrating 
in The First Longitudinal Mode 


48 



Using the one -force -per -half-wavelength concept illustrated 
in Figure 24, we can now proceed with the derivation of a value for the 
mode factor C m which relates displacement and force during a reso- 
nant condition. The following assumptions will be made: 

(a) The shape of the longitudinal bellows modes ie inter- 
mediate between a linear and a sinusoidal form. This 
is a very good assumption as comparison of the exact 
mode shapes with this approximation has shown. 

(b) The maximum deflection for a given mode is equal to 
the static deflection times Q, the dynamic amplification 
factor. 

(c) The maximum stress point for a given vibration mode , t 
occurs at the convolution with the maximum relative 
displacement* this is the end convolution for each half-wave 
mode shape, in general. 

The first assumption noted above implies that the mode shape is of the 
form illustrated in Figure 25, or over the first quarter wavelength 

*o ( ( N \ /Nrry\\ 

x = 2 }u ;,y + si n t )) < 17 > 

where x denotes the axial absolute displacement of a given point along 
the bellows defined by the axial position coordinate y; N is tin mode 
number. 

The second assumption given above, when considered in light of 
the one -force -per -half -wavelength concept, implies that the maximum 
absolute displacement x Q (see Figure 25) is 

FQ 

x 0 = 4Nk A (l 8 ) 

where F is the amplitude of the vortex shedding force applied as was 
illustrated in Figure 24. Again, Q is the dynamic amplification factor 
and K a is the overall bellows spring rate. 

From the third assumption given above, the point of maximum 
relative displacement x mr for a given bellows, and a given mode of 
vibration, occurs at a point where y = l /2N C . Therefore, from 





Equation (17), 


x o ( N , (v N\1 
x mr = 2 |n c t sm \2 N c / ) 

Combining Equations (18) and (19) yields 

FQ f N . Ar _N_\) 

x mr = 8 NK a |n c 8in \2 N c /) 

and comparison with Equation (16) shows that the "mode factor C m " 
defined in Equation (16), or 


now becomes 


x 


C m FQ 




( 21 ) 
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r = - — 
m 8N 

ITT” + sin 

U n J, 


( 22 ) 


and refers to the maximum relative displacement point (maximum stress 
point) along the bellows. 

At this stage, we are about in a position to estimate bellows flow- 
induced vibration amplitudes using Equations (15), (21) and (22), plus the 
data from Figure 22. We still need, however, vakes of Q, the amplifica- 
tion factor. A. comprehensive discussion of bellows damping and the 
resultant Q values will be delayed until later section ot this report; for 
the present we will assume the damping values are known. 

Test Data for Real Bellows 

The analysis method described in the foregoing has been used to 
correlate data from a number of flow-induced vibration tests of bellows. 
These tests were conducted in the following manner: 

(a) Strain gages were Installed on the tips of the end convolutions 
of all test bellows. On all bellows, one strain gage was 
oriented to monitor longitudinal strain and, on a few bellows, 
A gage to monitor circumferential strain was also installed. 
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(b) Each test bellows was subjected to a static load test, 
during which both deflection and stiain were monitored. 

This allowed force-deflection and force-strain curves 
to be plotted for each bellows. 

(c) Ea :h bellows type was subjected to a mechanical vibration 
test (to be discussed later) and the damping was determined. 

(d) Each test bellows was installed in the water-flow system, 
and flow-induced stress levels were monitored as a function 
of flow velocity. 

(e) This data was finally correlated with the analysis described 
earlier, 

The following discussion describes these various steps, and the results, 
in more detail. 

A typical strain gage installation on a bellows is shown in Figure 26. 
For most cases 1 /32-inch gages were used, while some 1 /64-inch gages 
were installed on those bellows with small width convolutions. Actual 
gage installation was straightforward, however, particular care was needed 
in lead installation to prevent fatigue. It was found that very small lead 
wire worked best (say about 0.005-0.010 inch) and that the best procedure 
was to glue the wirf>s down along the sides of the end convolution with 
rubber cement. The amount of cement used was kept to a minimum to 
ensure that the damping of the bellows was not significantly increased. 

Following strain gage installation, each bellows was tested to 
obtain force-deflection and strain-deflection curves. The reasons for 
doing this were twofold; first, it was desirable to have experimental spring 
rate values to allow more accurate mode frequency calculations and, second, 
by having a strain-deflection "calibration" along with the spring rate value, 
the flow-induced strain amplitudes could be directly interpreted as vortex 
shedding forces after suitable data reduction to account for the dynamic 
amplification, and mode factor. By this procedure, errors inherent in 
trying to analytically relate strain (or stress), and vibration amplitudes 
and forces, were eliminated. This allowed a more realistic judgment of 
the validity of the bellows flow-induced vibration analysis presented earlier 
in this report. Some typical force-deflection and strain-deflection data 
are shown in Figures 27 and 28*. 


* Appendix A gives dimensional data for all belloww used throughout the 
study; refer to this appendix for data on numbered bellows in all figures. 
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Deflection ( inches ) 






Strain ( it in. /in. 








The next step in the bellows test procedure was to subject one 
bellows of each particular type to mechanical vibration and thereby obtain 
a forced response plot for various conditions (various internal media, 
various internal pressures, different peak strain levels, etc.}. This 
mechanical vibration test procedure will be discussed in much mere 
detail in a later section. The important point here is to realize that 
the result of these tests was a tabulation of G values for each bellows 
for the various test conditions. 

The final step in testing each bellows was to mount it in the water 
flow loop and monitor vibratory strain levels as a function of internal 
fluid velocity. Figures 29, 30 and 31 illustrate typical flow-induced 
strain data for three different bellows types. In Figure 29, data for 
three supposedly identical test bellows is given, however, variations 
in vibration characteristics are readily evident. The mafn reason for 
this difference in the vibration characteristics is that the force -deflection 
and strain-deflection calibrations of these supposedly identical bellows 
were quite dissimilar, also. As will be shown, however, correlation 
of the apparently dissimilar data in Figure 29 for "identical* 1 bellows 
using the flow-induced vibration analysis discussed earlier gives very 
consistent results. 


This correlation for all data was carried out in the following 
manner: 

(a) The peak flow-induced strain data, as exemplified by 
the data in Figures 29, 30 and 31, was converted to 
equivalent maximum relative convolution displacement 
(Xm r ) values. The equation giving this conversion is 

x mr = ^ (dS/cLs)" 1 (23) 


\ 


where S is the measured peak strain and dS/dx is the 
slope of the overall bellows strain-deflection curve 
(see Figure 28). 


(b) Next, these maximum relative displacement values were 
used to calculate the quantity FQ from Equation (20), or 


8NK 


FQ = 


A 


mr 


Here, F is the peak vortex shedding force and Q is the 
dynamic amplification factor. 
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Dynamic Longitudinal Strain ( 















(c) From the results of the bellows damping study (to be 
discussed in a later section) the value of Q was known 
from the conditions of the test, therefore, the vortex 
force F was readily found from the results of (b) above, 

(d) Finally, the value of vortex shedding force F was used 
to calculate an effective vortex force coefficient Cp from 
Equation (15), or 


° F = (l/2pV 2 ) (A p ) 

where 


(25) 


Ap = ir(D^h + h 2 ) 

Figures 32 and 33 show the result of this data reduction process 
as applied to the flow-induced strain data of Figures; 29 and 31. The 
independent variable in Figures 32 and 33 is the Strc>uhal number based 
on convolution tip width. Note how much better the data fcr the supposedly 
identical bellows is correlated when reduced to the form shown in Figure 
32 than it was in Figure 29. 

The flow-induced strain data from all test bellows was reduced to 
the form shown in Figures 32 and 33. Because of differences in the con- 
volution geometry, the peak Cp values varied somewhat from item to 
item. In general, the greater the ratio of pitch to convolution tip width 
( X to ), the smaller the peak Cp value. This trend was also seen in the 
results of the single convolution test model discussed earlier, see Figure 
22. Unfortunately, the bellows tests discussed above did nc \ provide the 
range of convolution geometry required to adequately determine Cp as a 
function of geometry. Therefore, a special test procedure was developed 
to allow a better study of geometry effect. 

Bellows Geometry Effects 

Several tests have been performed in which the effect of bellows 
convolution geometry on the vortex force was determined. This was 
accomplished in a rather unique manner by step-wise stretching a^d com- 
pressing the test bellows so as to alter the convolution pitch. The advantag 
of this method are at least twofold; first, it allows geometry to be varied 
rather independent of other property variations, and second, it results in 
a hardware savings since more data can be obtained from a single test item 


Vortex Force Coefficient Cp 








Figure 34 shows a series of photographs of a sectioned bellows 
mounted in a stretching fixture and positioned in three different pitch 
configurations. Figure 35 gives flow-induced strain data taken from a 
bellows identical to the one shown in Figure 34. The data represents the 
observed peak flow-induced strain for each pitch setting. Figure 35 also 
shows this same data reduced to the vortex force coefficient form. Notice 
that the force coefficient generally decreases as the pitch increased; this 
trend is in agreement with the single convolution test model results (see 
Section II. 6) and otheV bellows tests. Also, however, Figure 35 shows an 
obvious trend of "optimum" pitch values which produce minimum local C-p 
values. This trend has been observed in several bellows tests and can 
probably be explained by vortex reinforcement as noted in Figure 9. These 
"optimum" pitch values can be correlated with the parameter k/cr and are 
most pronounced for the first few modes of a free bellows. The higher modes 
of a free bellows and flex hose modes don't seem to have these "optimum" 
pitch values; they apparently operate in a "vortex reinforcement condition" 
for all \/a values by a natural, and slight, "tuning" of the shedding frequency. 

Summary of Force Coefficient Data 

In the foregoing sections, vortex shedding force coefficient data 
has been presented as derived from essentially three different kinds of 
tests. First, the results of the idealized single convolution model tests 
were given in Section 11.6. These tests yielded force coefficient data for 
two different convolution configurations; the two configurations had the 
same convolute height h and the same tip width, but different pitches. 

Second, complete force coefficient data was obtained from tests of a 
number of "calibrated" bellows having various convolution geometries. 
Finally, a limited amount of test data was obtained by a bellows stretch- 
compression pitch changing scheme. Comparison of all of this force 
coefficient data shows consistent results which gives credence to the 
bellows flow-induced vibration model developed in the previous sections. 

Figure 36 shows two different curves representative of all available 
force coefficient data which we have obtained from the various sources. 

One curve is for use with the first few modes of a free bellows, and reflects 
the apparent vortex reinforcement trend discussed above. The other curve 
is for t!" i higher longitudinal free bellows modes, flex hose modes and the 
convolute local bending mode. 

This data represents, then, our state-of-the-art knowledge of vortex 
shedding force coefficients for bellows. It must be emphasized that this 
data can only be considered strictly valid for liquid flows through bellows of 
a size not too different from those tested to date (1.5 and 2.0 inch). We do 
not anticipate that these coefficients will vary much, if any, for other size 
ranges, however, the user of this data should be aware of a potential limi- 
tation. Larger bellows will be tested in the very near future to extend our 
knowledge of vortex force coefficients. 
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Figure 36. Summary Of Bellows Vortex Force Coefficient. Experimental Data 




II. 8 Bello ws Damping 


An extremely important factor, which provides a limiting condition 
i 1 the amplitude of the vibrations resulting from bellows flow excitation, 
is the damping of the coupled bellows -fluid system. The vortex force 
and mechanical model concept uiscussed earlier and illustrated in Figure Z4 
includes no damping so that, ideally, the vibration amplitudes would be 
infinite. Some damping is present, however, in all real structures and 
it serves to limit the vibration amplitudes to finite values. 


By way of review, consider the elementary linear second-order 
system shown in Figure 37. The forced response amplitude of this system, 
subject to a harmonic excitation F q sin c t, is of the form 


x(w) 




1/2 


(26) 


where 

w Q = (k/m) 1 ^, 


t, =b/2(k m )" 1/2 


At resonance (the frequency for which x(w ) is a maximum, the 

amplitude may be written as 
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In terms of the values of the mechanical elements of the system in 
Figure 37, Q may also be written as 
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Figure 38 illustrates a typical response curve for the system 
shown in Figure 37 where the amplitude of the force F q has been 
held constant, and the frequency of excitation varied* As shown, the 
value of Q physically represents the ratio of the vibration amplitude 
at resonance to the vibration amplitude in the limit of zero frequency. 
Also* the value of Q may be expressed in the form 

Q - Jl (3 0) 

Af 


where f r is the system resonant frequency and Af is the bandwidth at 
the halfpower point (0.707 times the resonant amplitude); see Figure 38. 

Because the value of Q is so important in the prediction of 
bellows flow-induced vibration amplitudes, and even for proper data 
reduction as discussed earlier, a large number of bellows mechanical 
vibration tests were conducted from which Q values could be obtained. 
Figure 39 shows the test s^t-up used for the mechanical vibration tests. 
Each bellows was mounted in a special fixture which rigidly connected 
the end flanges. The reason for this was to ensure that the vibration modes 
excited would be the same ?,s those observed in the flow-excitation tests. 


Figures 40 and 41 show Q values obtained from the forced vibration 
tests. In general, the following trends are noted (see also Figure 42): 

(a) The bellows damping is generally most sensitive to peak 
strain amplitudes, with the damping increasing as strain 
increases . 

(b) Single -ply bellows of different configurations, tested under 
similar conditions (same strain level, same internal media, 
etc.) do not necessarily have the same Q values. 

(c) For one configuration of bellows, increasing the number of 
plies from one to two significantly increases the damping 
(reduces Q); however, going to three plies adds little more 
damping, if any, over the two-ply configuration. 

(d) For bellows with a liquid as the contained media, damping 
may greatly increase at low absolute pressure's. 
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Figure 38. Example Bellows Forced Response Data For Constant Level Acceleration Input 
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The fact that damping increases as the strain level increases 
means that the damping is a nonlinear function of amplitude and, 
therefore, the linear model in Figure 37 is not necessarily valid. 

Valid Q values can, however, still be assigned if it is realized that 
they depend on strain level . The fact that the bellows damping is non- 
linear is really no surprise, since structural damping is often of this 
nature (see Lazan's discussion of hysteretic damping in Reference 11). 

The second trend noted above is that bellows with different con- 
figurations do not necessarily have the same Q values for similar test 
conditions. This is undoubtedly caused by a lack of a proper correlation 
model for the bellows damping. In other words, if we could develop an 
analytical model which allows prediction of damping in terms of convo- 
lution configuration, strain level, etc., then the damping data should 
’’correlate". Such a model has not yet been developed, but efforts along 
these lines will be attempted. For the present, we have no recourse but 
to use the data in Figures 40 and 41 to define conservative Q values for 
design and analysis purposes. 

Table V and Figure 43 summarize what we consider to be a best 
estimate of conservative Q values for various situations. The user of 
this data is advised to always tend toward the higher Q values where a 
choice must be made. 



II. 9 Other Influences 

In practice, it has been found that several factors, in addition to 
those already discussed, can have an influence on the magnitude and 
character of bellows flow-induced vibrations; these factors are, (a) 
presence of an upstream elbow, (b) angulation of the bellows, (c) non-rigidly 
attached piping, and (d) acoustic resonances of the duct system. The first 
three influences are discussed briefly in this section and the subject of 
"acoustic resonance effect" is treated in the following section. 

Upstream Elbow Effect 

Because it had been observed that some bellows installed immediately 
downstream of a duct bend were seemingly more susceptible to flow-induced 
failure problems than those in straight sections, a small investigation was 
undertaken to study this. One of the calibrated test bellows was placed in 
the flow system, as shown in Figure 44, and located immediately down- 
stream of a short-radius 90-degree elbow. Flow-induced strain ampli- 
tudes were then monitored as a function of the mean internal fluid velocity. 
The test was repeated with the same bellows in a straight run with all 
other conditions kept the same. 
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Specific Spring 

Number 

Internal Media 

Curve 

Rate (see Note 1) 

Plies 

(see Note 2) 

No. 

all ranges 

i 

low pressure gases 

1 

over 2000 lb/in^ 

i 

high pressure gaces, light liquids 

1 

over 2000 

i 

water, dense liquids 

2 

under 2000 

i 

high pressure gases, light liquids 

2 

under 2000 

X 

water, dense liquids 

3 

over 3000 

2 

all 

3 

2000-3000 

2 

all pressure gases 

4 

under 2000 

2 

all pressure gases 

5 

2000-3000 

2 

all liquid s 

5 

under 2000 

2 

all liquids 

j 

6 

over 3000 s 

3 

all 

4 

2000-3000 

3 

all 

5 

under 2000 

3 

all pressure gases 

5 • 

under 2000 

3 

all liquids 

6 


Use of Table - To use table, first calculate bellows specific spring rate, 
then look up application curve number corresponding to 
this specific spring rate, number of plies, and internal 
media. 

Note 1: The specific spring rate is here defined as 

S.S.R. = k a! 1 c - 

or is the spring rate per convolute, per ply, per unit of 
diameter. 

Note 2: Low precsure gases will be defined here t-s being those 

gases below 150 psia. Light liquids will be defined as 
having a density, relative to water, of less than 0.2. 


TABLE V 

Applications Information for Use with Q Values 
Data in Figure 43 
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Figure 45 shows the results of these tests. Notice that two things 
have changed about the response of the bellows with the elbow present, 
compared with the response with no elbow. First, the mean fluid velocity 
range over which flow excitation occurs has been shifted down by about 
45 percent. Second, note that the strain level for a. given mode of 
vibration has significantly increased, or that the strain amplitude at a 
given fluid velocity has increased. 

The explanation for the difference in bellows vibration with, versus 
without, an upstream elbow is that the velocity profile is skewed or dis- 
torted because of the presence of the elbow (see Figure 45). This causes 
an excessive dynamic fluid loading on the convolutions downstream of the 
outer bend region and possibly a reduced loading downstream of the inner 
region; the net effect is an increased effective value of the fluid velocity 
along the surface of the bellows, causing the vortex shedding to lock-in 
on the bellows vibration at a lower -than-normal mean fluid velocity. 

While further tests must be performed to better quantify the 
influence of an upstream elbow on bellows vibration levels, we can say 
for the present that an increase in stress level by a factor of about Z is 
possible for a given mean velocity. Elbows with more generous radii 
of curvature, or located further upstream, would most likely not cause 
this great an increase in flow-induced stress levels, 

Non-Rigid Ducting 

As might be expected, if the ducting attached to the ends of a 
bellows is not rigid then the character of the flow-induced vibration can 
be changed compared with the ideal case of perfectly rigid attachments. 
There are apparently two ways that the bellows vibrations can be affected; 
first, the modal frequencies can be changed (generally lowered), and 
second, additional damping can be "added” or energy dissipated through 
the attachments. 

For most cases, the modal frequencies of a bellows will not be 
significantly changed because of flexible attachments. The reason for this 
is that the bellows is generally so much more flexible than the attached 
ducting, particularly in the longitudinal or axial direction. Apparently, 
however, significant damping can be introduced by allowing duct flexi- 
bility as evidenced by some of our experimental observations. While we 
can’t at this time report any quantitative results, the reader should take 
note of these observations. Whenever a bellows is undergoing a qualifi- 
cation test, our observations indicate that particular care should be 
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Figure 45. Belows Fiow- induced Strain With And Without Upstream Eibow 
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exercised so that the bellows attachments are at least as rigid as in 
the actual installation. This should help ensure that the flow-induced 
vibration levels attained in the qualification test are at least as great 
as the actual installation. 

Angulation of Bellows 

At the present time, we have no quantitative information to 
report concerning the effect of angulation on bellows flow-induced 
vibrations. It would seem, based on the results of the upstream elbow 
tests discussed above, that angulation might lead to an increase in 
vibration levels because of the excess fluid loading on the outer bend 
region. Also, however, angulating a bellows might cause a change in 
the mode frequencies and damping. At this time, effect of angulation 
must be considered unknown and we caution the reader to be conserva- 
tive in dealing with this situation. 

II. 10 Acoustic Resonance 


In a previous section entitled "Fluid Leading Effects, 11 we 
discussed how acoustic loading can increase the number of degrees of 
freedom of a bellows system by producing coupled bellows-acoustic modes. 
When this type of coupling occurs, it has been generally observed that 
very significant pressure fluctuations can be propagated throughout the 
fluid system. Appendix D discusses a mechanism of noise production 
as a result of a vibrating bellows exciting the longitudinal acoustic duct 
modes. In addition, however, acoustic resonances have been observed 
to cause a significant increase in the flow-induced stress levels. This 
second effect is apparently a result of the acoustic resonance flow and/or 
pressure fluctuations coupling with the vortex shedding process to pro- 
duce a force amplification, as will now be discussed. 

In the preliminary air flow testing of one free bellows, it was 
noted that the specimen showed virtually no response up to a rather high 
velocity where, quite suddenly, the response became very significant. 

This sudden onset of response was accompanied by very significant 
acoustic radiation from the entire flow loop, indicating that an a.coustic 
resonance existed in the facility. A calculation of the stress indicator 
(discussed in the next section) value at the air velocity where this 
phenomena started, gave a value of only a few thousand psi. It was felt 
that the vibration observed was much more severe than it should have 
been, and this was attributed to an enhancement of the vortex shedding 
by the acoustic resonance. At this point it was believed to be a longi- 
tudinal resonance which caused this effect. 
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Following this preliminary air flow test, two different bellows 
were extensively investigated in order to find out more about the 
influence of acoustic resonance on the bellows response. One of these 
bellows was tested using two different gases; air and Freon 12, These 
two gases were selected because of the great difference in speed of 
sound. Under ambient conditions, air has an isentropic speed of sound 
approximately equal to 1100 feet per second, while the isentropic speed 
of sound for Freon 12 is about 485 feet per second. Also, the Freon has 
a greater density (by a factor of about 5) under the, same conditions. To 
be able to assess separately the effect of density and speed of sound, 
flow tests were conducted at different mean pressure levels with each 
gas. Therefore, by comparing the test results with one gas at different 
pressures, the effect of density was determined. Figure 46 shows, for 
example, acceleration level at the predominant response frequency, 
(measured at the bellows flange) as a function of mean air velocity for 
two different mean pressure levels. Qualitatively, the response is the 
same for the two cases, and shows no significant vibration up to a critical 
flow velocity, where, suddenly the bellows subsequently responds well in 
various flow-induced modes of vibration as the velocity is further increased. 
This critical velocity is the same for both mean pressure levels, as may be 
seen in Figure 46. This onset of bellows vibration is accompanied by a 
great increase in the acoustic noise level at the same predominant fre- 
quencies that exist in the flange acceleration signal. 

The results of the tests conducted with the Freon were qualitatively 
the same as with the air; the critical flow velocity at which the bellows 
suddenly began to respond was, however, much lower. In fact, the 
critical flow velocity with the Freon was about half that experienced with 
the air. This is just what we expected might happen, assuming the 
responsible mechanism to be an acoustic resonance, since reducing 
the speed of sound would cause a reduction in the frequency (hence, flow 
velocity) required to give a resonant condition. 

To further verify the acoustic resonance idea, the acoustic pres- 
sure field inside the bellows was mapped out, with the aid of a small 
microphone; the result is illustrated in Figure 47. It was determined 
that we were getting excitation of coupled radial-axial acoustic modes, 
with the maximum sound pressure level, in the radial direction, occurring 
in the vicinity of the bellows internal diameter, while the sound pressure 
level maximum, in the axial direction, is half way along the convoluted 
length of the bellows. There is an acoustic pressure minimum all along 


82 




Fluid Velocity - ft/ sec 2209 
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Figure 47. illustration Of Acoustic Pressure Level Variation 
With Radial And Axial Position In Bellows 
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the geometric bellows axis. The acoustic field described above is 
being excited by the vortex shedding from the bellows convolutions. In 
turn, the acoustic resonance apparently enhances the vortex shedding, 
causing "greater-than-expected 11 response levels. 

The reason for the sudden occurrence of this phenomena at some 
critical flow velocity is as follows: At low fluid velocities, hence low 
vortex frequencies, the acoustic wavelengths are too long to "fit" in 
the bellows in a radial direction; therefore, the radial acoustic waves 
which are generated experience pure attenuation. When a high enough 
shedding frequency is achieved (hence, a high enough fluid velocity), 
the wavelengths are short enough to "fit" radially, and, therefore, cause a 
radial resonance condition and couple with an axial acoustic mode. 

It is very important to note that this radial acoustic resonance 
was observed even when the bellows was completely restrained to prevent 
any possible convolution vibrations (this was done by covering the bellows 
on the outside with a thick coat of epoxy). Under these conditions, the 
acoustic resonance was not as strong as with the bellows free to vibrate. 
Even so, accelerometer data monitored at the bellows flange showed 
strong acceleration levels at the acoustic resonant frequency, plus 
numerous harmonics. This implies that the use of an accelerometer 
for bellows vibration monitoring where an acoustic resonance occurs 
can be very misleading. A comparison was made during one test of 
the spectral content of convolute strain and flange acceleration, and 
the results were grossly different. The reader should, therefore, be 
wary of accelerometer data which is taken where an acoustic resonance 
occurs, and which shows strong response at multiples of the vortex 
shedding frequency. 

In conclusion, the effect of a radial acoustic resonance in bellows 
leads to greatly increased flow-induced stress levels and sharp increases 
in the sound level in the duct system a v the vortex shedding frequency, 
and at multiples of this frequency. A rough judgment can be made as to 
the likelihood of a radial, acoustic resonance occurring by comparing the 
calculated acoustic wavelength with the outside diameter of the bellows. 
When this calculated wavelength becomes small enough to about equal the 
bellows outside diameter, then acoustic resonance may occur, and any 
bellows structural modes corresponding to this frequency may be strongly 
excited. 



II. 11 bellows Flow -Induced Stress 


Best Model 

The present n best model 11 for predicting bellows flow-induced 
stress levels involves calculating the maximum relative convolute 
displacement by the method discussed previously in this chapter, and 
then calculating a corresponding stress level from an existing stress- 
deflection theory. The steps in thus calculating flow-induced stress are 
summarized below: 


Step A - Calculate the maximum convolute relative displacement 
of one -half of a convolute from the expression 


x = 


Cm F Q 

Ka 


(31) 


where 
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= J_f 

8N ( 
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iC + Sln V2 nT)1 


(32) 


and 


F = C F tt (Dih + h 2 ) (l/2pV 2 ) (33) 

The values of Cp and Q may be found from Figures 
36 and 43, respectively. The mode number N may 
be determined as a function of fluid velocity by the 
method outlined in Section II.4(see Tables III and IV). 

Step B - Calculate the stress corresponding to the maximum 
convolute displacement x by an existing method, 
such as outlined in Appendix B. 

Step C - If an upstream elbow is present for the subject bellows 
installation, then some upward adjustment of the calcu- 
lated stress is necessary. For a severe bend immediately 
upstream of a bellows, multiply the calculated stress by 
2.0 to correct for the elbow. For "less severe" cases, 
multiply by 1.5 to correct for the upstream elbow. 


At this point it is well to point out the currently known limitations 
and potential limitations of this model, including estimated possible errors. 
The model, as discussed in this chapter, has evolved from a theory which 
is supported by some empirically- obtained factors; these are the vortex 
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force coefficient Cp and the amplification factor Q, Values of C F 
and Q were obtained from tests of bellbws, and a special single convolute 
model, having internal diameters of 1 . 5 inches and 2,0 inches. Bellows 
with a range of geometries were tested in these two sizes. It is evident, 
therefore, that the model as it stands has only been verified for bellows 
of this general size. A potential limitation, then, is in the use of the Cp 
and Q values reported herein for much larger (or smaller) bellows. Also, 

Cp data was only obtained with internal water flow. However, we are 
confident that the Cp values thus obtained can be used for gas flow through 
this general size bellows i f no acoustic resonance occurs . 

Probably the greatest potential source of error in calculating flow- 
induced stress for a given bellows is the value picked for the amplification 
factor Q from Figure 43. As discussed previously, we have not yet 
developed what we believed is a satisfactory damping correlation method 
which includes the effect of all factors illustrated in Figure 42, The data 
given in Figure 43 is felt to be conservative, or give !! liigh n Q values. 

Other sources of error are the occurrence of acoustic resonances 
in the duct system and non-rigid M end conditions 9 * 11 . Acoustical resonances 
tend to amplify the vortex shedding forces, while non-rigid attached ducts 
tend to reduce vibration amplitudes because of increased damping. 

The sources of error and their potential effect on the calculation 
procedure are summarized below: 

(a) Vortex Force Coefficient Cp - The Cp values given in 

Figure 36 are conservative in that they represent the upper 
limit of the experimental data. Our observations show that 
these Cp values are good for bellows vibrating vigorously 
so that there is M good coupling n between the vibrating bellows 
and the vortex shedding. For bellows which do not vibrate 
vigorously (have low stress values) these Cp values will tend 
to be high; but this will result in conservative answers. 

(b) Dynamic Amplification (Q) Values - The values of Q given 

in Figure 43 are conservative, or at the upper limit of experi- 
mental observations. For a very "stiff" bellows, these values 
should be quite good, but for a "low spring rate” bellows they 
will tend to be high but give conservative answers. 

( 9 ) Acoustical Resonance (Radial and Helmholtz) - An acoustical 

/ resonance of the type discussed in Section II. 10 will tend to 

make the calculated stress levels very conservative . As a 
rule of thumb, we suggest that the acoustical resonance may 
cause a stress amplification of the order of 5. 
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(d) Upstream Elbow - The presence of an upstream elbow 
causes excessive flow-induced stresses for a given fluid 
velocity. We don't have much quantitative data on this 
effect but suggest correcting the calculated stresses 
upward by as much as 2.0 to account for the elbow. 

(e) Heat Transfer Effect - A potentially great influence on 
bellows flow-induced stress can be caused by heat transfer 
induced changes of phase or fluid state within the bellows. 
This effect is currently under investigation. 

Improved Stress Indicator 

In Quarterly Report No. 5 (16) , we developed the preliminary 

form of a quick method for estimating the dynamic flow-induced stress 
level for a given bellows installation. This method consisted of calculating 
a "Stress Indicator" which was a number roughly proportional to stress and 
apparently valid for judging the relative severeity of bellows vibrations. 
We are still convinced that the "Stress Indicator" approach is useful and 
valid. The purpose of this section is to present an "Improved Stress 
Indicator" which is based on the "best model" stress calculation procedure 
discussed above. 

The basic equations which describe flow-induced stress, as given 
in Sections II. 5 to II. 7, are 

F = C F Ap(l/2 pV 2 ) = C F A p p d (34) 
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C m K2 

Ka 


and 


Stress = 


C s Etx 

h 2 


(35) 


(36) 


The objective of the "Stress Indicator" approach is to combine these 
expressions in a simple form so as to give a single expression for stress 
which contains only readily-known bellows dimensional data, or parameters, 
and flow variables. Therefore, we will assume the following: 



(a) 


A can be approximated by 

IT 


(b) 

(c) 


A p * TrD m h 

Cj. is known from Figure 36 
C m can be approximated by 


C » — — 
m 4N. 


(d) 


can be approximated by 


K A * D m E <V N c> 

(e) Q can be found from Figure 43. 

(f) The factor C g is approximately constant for all tallows. 

Based on these assumptions, Equations (34), (35), and (36) maybe combined 
to give 


/C F QWh\2 

Stress = irC s ^ jj J j (1/2 pV 2 ) (37) 


Since C g has been assumed constant, we can extract a quantity from 
Equation (37) which is approximately proportional to flow-induced stress 
for all bellows; this quantity is 


Stress Indicator 
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(38) 


This expression for Stress Indicator is slightly di f ferent from that given 
in Reference (16) and is closer to the real situation. To perform calcu- 
lations with Equation (38), it is first necessary to find Cp and Q from 
Figures 36 and 43, respectively. 


Bellows Fatigue Life 


To supplement the Stress Indicator given in Equation (38), we have 
prepared a preliminary plot (Figure 48) for predicting bellows fatigue life 
which gives the number of cycles to failure as a function of Stress Indicator. 
The actual time to failure can be calculated by dividing the number of cycles 
by the appropriate modal frequency (see Section II. 4, Tables III and IV). 
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Notes : 


- Obtain Cp from Figure 36 i j 

- Obtain Q from Figure 43 j j 

- Use Cp 3 2,0 for severe upstream elbow 

- Use Cp = 1.0 for no upstream elbow 

- This data is unconservative where acoustic resonance occurs 
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Figure 48. Preliminary Bellows Fatigue Life Data 


The data in Figure 48 represents a collection of all failure cases 
which we have observed in the laboratory, plus some data for the F-l 
engine ASI line failures, and for S-IVB J-Z engine LOX feedline bellows 
(with appropriate accounting of the upstream* elbow). This fatigue data 
must be considered preliminary because of the small number of sample 
points used in its compilation, but should be good enough, at the very 
least, to help judge relative severity of various bellows installations. 
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III. BELLOWS LINER DESIGN AND EXTERNAL 
DAMPING DEVICES 


II1.1 Cone Liners for Reduction of Vibrations 


A limited study of cone-shaped liners was made with the specific 
purpose of determining how much of the bellows must be covered by the 
liner to effectively stop flow-induced vibrations. The liners being con- 
sidered here are simple, truncated, hollow cones formed from thin sheet 
metal which are mounted at the entrance of a bellows and extend into the 
bellows itself. The cone shape allows angular movement of the bellows but 
also produces a constriction, and a resulting pressure drop, which is 
discussed in the next section. 

An initial aspect of the liner study was to determine the nature 
of the velocity profiles downstream of the liner; this, it was hoped, would 
provide seme insight into how far the liner should extend into the bellows 
to prevent vibration. Pitot tube traverses were run downstream of con- 
vergent inserts mounted in a 2-inch ID duct to determine the shape of a 
fluid jet as it emerges from the end of a cone liner. The resultant 
velocity profiles are shown in Figure 49 where the velocity downstream 
has been normalized by the velocity at the center of the cone exit. 

The influence of the cone angle was not very pronounced, but there 
was a vena -contr acta present with the 30° cone which was not present with 
the 10° or 20° cones. This was quite noticeable from the Pitot differential 
pressure which represents the square of velocity. However, the actual 
center velocity at the vena-contracta was only about 3% greater than at 
the cone outlet. Comparison of the flow profiles at station #6 with and 
without the inserts shows that the flow is noticeably affected as far as 5 or 
6 outlet diameters away from the insert. This would indicate that a con- 
siderable part of a bellows could be left uncovered when a cone-shaped 
liner is used. 

In order to determine the extent to which this information could be 
applied to an actual situation, two -1/2 -inch ID bellows of different convolution 
geometries (#110 and #112)'* were flow tested with water using a movable 
liner that could be positioned at any axial location within the bellows. Since 


See Appendix A for a list of dimensional data on all test bellows. 





the effect of cone angle on the flow pattern was rather small. only one 
liner configuration was tested. The liner had an outlet diameter of 
1-1/4", a length of 11/16", and aT0° included angle. The fluid velocity 
was varied from 0 to 7 0 ft/ sec based on bellows diameter and maximum 
dynamic strain was recorded for several liner positions. The results 
are shown in Figure 50. 

Even with the liner completely withdrawn from the bellows, the 
strain level was reduced from that of an unlined bellows; this result is 
in agreement with the velocity profile observations. However, all modal 
vibrations were not stopped until the liner covered all but one or two 
convolution roots. The modes which were present did not occur in order. 
In some cases the first mode was present at two different flow rates while 
the second mode was completely absent over the entire flow range. 

Operating the bellows at an angle did not significantly affect the 
dynamic strain readings for these particular bellows -liner arrangements. 
There probably are cases where the effects of angulation are significant. 

The effectiveness of a cone liner or the effect of angulation would 
be difficult to predict analytically due to the highly turbulent expansion 
which takes place at the liner exit. The tests which have been performed 
can serve as a guide, but more testing of different types of bellows is 
required before reliable empirical design is possible. At this time, how- 
ever, we offer the following design guidelines: 

(1) If the bellows failure problem without a liner is severe, 
then the liner should cover most if not all of the active 
convolutions. If pressure loss is a problem, the liner 
might be reduced j k length somewhat, possibly allowing 
some bellows vibration. Testing should be performed to 
verify the design; preferably, a ’’before and after" test 
should be conducted with one or more bellows. 

(Z) If the bellows vibrations are not too severe without the 

liner, then a very short liner extending over, say, about 
Oi.e-fourth of the convoluted length should offer substantial 
reduction in vibration levels. Again, the design should be 
verified with testing. 

(3) On any liner design, an important consideration is weight 
which means the liner thickness, as well as length, should 
be minimized. Any proposed design should be subjected to 
a buckling analysis to assure that it will not collapse as a 
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Figure 50, Maximum Dynamic Strain For Fluid Velocities From 

0 -70 ft/sec (Water - 75°F ) 
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result of the differential pressure. Also, the designer should 
be aware of potential "splitting” problems caused by turbu- 
lent flow excitation of tne liner. High dynamic pressure 
situations should be approached with care. 

Ill . Z Pressure Loss of Cone Liners 

An important point for consideration in regard to a proposed 
liner design is that the pressure loss be acceptable. It has been demon- 
strated that a liner can, in fact, increase the pressure drop compared 
with the bare bellows. The purpose of the following discussion is to 
demonstrate a simple procedure for estimating liner pressure drop 
characteristics . 

Usually, a typical liner can be considered, from a fluid mechanics 
point of view, to consist of a contracting section followed by an expansion 
or, in some cases, an expansion followed by a contraction or, possibly 
even a series of expansions and contractions depending on the specific 
configuration. Assuming this to be true, then using information available 
in the literature for losses resulting from expansions and contractions, 
the liner overall pressure loss can be estimated. Figure 51 gives all of 
the information needed for estimating these losses. Note that for either 
a contraction or expansion, the pressure drop is characterized by 

A Kjl „ 2 - , ^ ,2 n 2 

Ap = 2 V 1 [1 -Pi/D 2 ) ] 

= if V 2 2 [(D 2 /D 1 ) 2 -1] 2 

where K is a loss factor defined in terms of the geometry of the duct as 
shown in Figure 51. 

This prediction procedure has been checked experirnehtally for the 
liner illustrated in Figure 52 which is a conical configuration typical of a 
currently used design. 

For a zero deflection angle, the pressure loss for this case will 
be the collective result of, first, the contraction in the cone section and, 
second, the uncontrolled expansion when the fluid exits the cone. The 
pressure drop will therefore be given by 

Ap = i (K C + K E )[(D 2 /D 1 ) 2 -!] 2 V 2 2 
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Calculations Of Pressure Drop 

For Either Contractions Or Expansions , The Pressure Drop Is 
Ap a Vf [l-(D,/D 2 ) 2 j ;j ~g~vt [( / D| ) 2 - lj 

Figure 51. Flow Loss Factors For Contractions And Expansions 
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From Figure 51, the value of K c for this cone geometry is about 0.35, 
The value of Kg for an uncontrolled expansion is 1.0, thus 


Ap =f (1.35) [(1.22) 2 3 -1] 2 V 2 2 = (0.324) (1/2 pV 2 2 ) 


Figure 52 shows pressure drop as a function of mean flow velocity predicted 
by the above equation, and also the pressure drop obtained experimentally 
for this same conical liner configuration. Figure 52 al^o shows the experi- 
mentally obtained pressure loss for the bare bellows (no liner) used for 
the tests. In this case, the liner is responsible for an increase in pressure 
loss. 

III. 3 External Damping Devices 

One possible means of reducing the stresses in a bellows resulting 
from internal flow is to add damping to the bellows structure. The attrac- 
tive feature of this method of reducing stress is that it allows the designer 
to possibly "fix" a bellows which has demonstrated some fatigue tendencies, 
but is otherwise satisfactory. Our work to date in this area has been 
strictly exploratory; several external damping ideas have been tried 
with two basic objectives, (1) to find out in a quantitative fashion how 
much the damping can be increased under ideal conditions, and (2) to 
compare the relative increase in damping for several devices or ideas. 
Three types of external damping ideas have been tried thus far; these are: 

(1) Filling the convolutions with some damping material, 
usually a viscoelastic material, such as RTV rubber, 

3M strip-caulk (a putty-like material), and a polystyrene 
foam. 

(2) Placing a metal spring or rubber "O" ring in each 
convolution. 

(3) Wrapping the exterior of the bellows with screen wire. 

Figure 53 shows a photograph of several bellows with various 
external damping devices added. Flow tests have been run on several of 
these "damped" bellows to explore very briefly their possibilities. Figures 
54 and 55 show some of the test results. In Figure 54, a comparison is 
shown between an undamped bellows and the same bellows with two types 
of dampers; the external screen and the neoprene "O" rings. In general, 
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both devices provided considerable damping; however, the n O M ring 
dampers allowed some resonance effect at about 60 ft/sec flow velocity, 
whereas the external screen did not. Figure 55 shows a comparison of 
an undamped bellows and two others with damping devices; one of the 
damped bellows had 3M strip -caulk packed in the convolutions and the 
other had coil springs wound around the outside of each convolution. 

The two damped bellows exhibited about the same vibration character- 
istics, as may be seen from Figure 55, and both reduced the vibration 
level by a considerable amount (about 2 orders of magnitude near the 
resonant points). 

In general, we believe that these exploratory tests have demon- 
strated that it is, indeed, possible to add considerable damping to a bellows 
by use of some external device. A word of caution must be added here, 
however, because our tests were only conducted on free bellows and for 
conditions where the first few modes could be excited. For a flexible 
hose or the higher modes of a free bellows, it is not known what kind of 
device will be required. Of course, some of the devices which were 
t^ied would obviously not v/ork for a flexible hose. 

Another point which the reader will undoubtedly note is that some 
of the damping devices which were tried are absolutely not practical for 
c.xe reason or another; for example, the external screen, as shown, 
would not allow a free bellows to flex; a variation of this idea will, how- 
ever. Also, neoprene, the strip-caulk, and the polystyrene foam will 
not work at cryogenic temperatures. As mentioned previously, one 
objective of these tests was to find out how much damping can be intro- 
duced under ideal conditions, and this is wfyy some impractical devices 
were tried. The important pomt of the results is that considerable , 
damping can be added; the problem now is to discover practical devices. 

III. 4 Summary of Bellows Vibration Suppression 


We have demonstrated in this chapter that both cone liners and 
external damping devices are effective in suppressing bellows flow- 
induced vibrations. Both types of devices have been used in the past by 
bellows manufacturers for this very purpose. Seme points are important 
to reemphasize about the use of these devices: 

(1) Cone liners can cause a significant increas e in the 
pressure drop across the bellows. 

(2) Cone liners can experience both buckling failures, 
because of the large differential pressure, and fatigue 
failures because of turbulence -induced vibrations. 
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(3) External damping devices can be as effective as a cone 
liner without causing an increase in pressure loss. 

(4) Care must be exercised in the application of external 
damping devices to ensure that the damping is actually 
increased. 


i 
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BELLOWS AND ELBOW PRESSURE LOSS 


IV. 1 Description of Pressure Loss Mechanism 

As described in Section II. 3, the general behavior of the fluid 
through a bellows involves a generation of vortices in the convolutions. 

In the absence of vibrations, the vortex motion inside each convolute 
is steady as was shown in Figure 6. With flow-induced vibrations, 
toroidal vortices are periodically shed out of the bellows convolutes as 
shown in Figure 9. 

Regardless of whether or not there is flow-induced vibration, the 
generation of vortices in the convolutes causes an extraordinarily high 
steady flow pressure loss compared with an equivalent length of smooth 
tubing. Figure 56 shows a typical plot of the friction factor f as 
compared with friction factors for flow through tubes. Recall that the 
friction factor is proportional to the pressure drop in a one-diameter 
length of tubing (or bellows), and is related to the pressure loss and 
fluid conditions by the equation 

Ap = (fL/D) (i/2 pV Z ) 

Note from Figure 56 that the bellows friction factor increases 
drastically as the Reynolds number is increased over the range from 
about 1C>4 to 10 6 . This trend is contradictory to the observed behavior 
of the friction factor for turbulent flow in a pipe, which shows that f 
decreases as Reynolds number increases. It is our present belief that 
the cause for this peculiar behavior can be explained by the presence of 
the vortices in the convolutes, and the fact that these vortices will trans- 
ition from a laminar to a turbulent condition at a higher Reynolds number 
than does the main stream flow. We are planning to explore this idea with 
a simple analysis in the near future. 

IV . 2 Compiled Bellows Pressure Loss Data and Some Test Results 

All pressure loss data for bellows and flex hose which was available 
in the open literature, or from unrestricted company reports, has been 
compiled; this data was presented in Quarterly Report No. 1 (Reference 14). 
This compiled data included two different correlation methods giving the 
friction factor f in terms of bellows geometry and Reynolds number; these 
two correlation methods are illustrated in Figures 57 and 58. Each of 
these correlations gives f, as a function of Reynolds number, in terms of 
a one-parameter family of curves. The parameter used in the Daniels 
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Figure 56. Comparison Of Friction Factors For Smooth Pipes And Flexible Hose 
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Figure 57. Bellows Friction Factor Correlation From Reference [ 12 ] 
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and Cleveland correlation (Figure 57 and Reference 12) is the ratio 
of convolute height to bellows internal diameter. The parameter used 
in the Riley correlation (Figure 58 and Reference 13) is the ratio of 
convolute pitch to internal diameter. Both correlations can only be 
valid at the same time for a series of bellows having convolutions with 
exactly similar geometry; i.e., the ratio of convolute height to pitch, 
etc., remains constant for all bellows of the series. Therefore, one 
or both cannot be completely valid for all bellows, in general. 

To get a better insight into bellows pressure loss as a function of 
geometry, and to attempt to clear up the discrepancy of correlation 
methods noted above, a series of experiments were conducted with some 
selected bellows, and with a special convoluted test model. The special 
test model, described in Quarterly Report No. 2 (Reference 15), had 
provision for changing the convolution geometry so that the effect of 
geometry on pressure loss could be determined in a controlled manner. 
Figure 59 shows some pressure drop data obtained with this model. 

Figure 60 shows friction factor data reduced from the pressure drop 
data of Figure 59. All of the data in Figures 59 and 60, except one case, 
were for simulated convolutions with identical depths (or heights). Com- 
parison of the data for this condition (curves 1, 2, 3 and 4) show a distinct 
variation in the friction factor. This data tends to invalidate the correla- 
tion given in Figure 57 since it shows a variation in f for constant values 
of h (the ratio of convolute depth to internal diameter). Similarly, a 
comparison of curves 1 and 2 in Figure 60 shows some variation in f for 
bellows having the same value of X/ (the ratio of pitch to internal 
diameter). This tends to invalidate the correlation given in Figure 58. 
Thus, we can conclude that neither of the correlations given in Figures 57 
and 58 is entirely or generally valid. Further detailed evaluation of these 
two correlation methods has been performed by comparison with more of 
oar experimental data and that of other investigators. We have concluded 
that: 

(1) The friction factor f is almost entirely independent 
of h/D^ (the ratio of convolute height to internal 
diameter); therefore, the correlation given by Daniels 
and Cleveland (Figure 57 and Reference 12) has no 
meaning for bellows of general geometry and is invalid. 

(2) The correlation given by Riley (Figure 58 and Reference 13) 
is roughly valid for all bellows. However, the, Riley -type 
correlation could be improved by making it a three - 
parameter correlation, with the three parameters being 

X/D^, \ /a and the fluid density. 
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Figure 59* Pressure Drop As A Function Of Flow Rate 
For Five Beliows Test Geometries 
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We recommend that the Riley correlation be used for pressure 
loss estimate purposes with the user aware of the fact that this method 
generally overpredicls pressure loss for liquids and underpredicts for 
gases. Also, the Riley correlation tends to overpredict for bellows 
with "open" convolution designs ( X./cr > 2.0), and underpredicts for 
bellows with "closed"convolution designs ( \/a < 2.0). 

In very general and rough terms, the trends of "good" and "bad" 
convolutions configurations are illustrated in Figure 61. It is interesting 
to note that the "bad" designs correspond to those geometries which exhibit 
the more severe flow** induced vibration tendencies. 

IV. 3 Low Pressure Loss Bellows Configurations 

As discussed in Chapter III, a conical liner can be effective in 
suppressing flow-induced vibrations but can lead to a significant increase 
in the pressure loss if the required contraction is too great. For bellows 
not required to angulate, a cone liner can reduce the pressure loss if it 
is property designed. Therefore, new ways were sought to reduce losses 
which might be useful for bellows required to angulate. 

Several specific configurations have been tested. A sketch of one 
configuration is shown in Figure 62. This design employs a rigid skirt, 
but not one of conventional design. Rather than neck down the skirt or 
liner section, the convoluted section is made slightly larger (on the order 
of 10%) and the skirt is kept the same diameter as the duct. The ratio of 
skirt diameter to bellows diameter is, of course, dictated by angulation 
requirements. Just downstream of the bellows a converging section 
guides the fluid into the downstream duct. Figure 63 shows pressure loss 
as a function of flow rate for this new design as well as for a conventional 
skirt design, and unskirted bellows and a smooth wall section of tubing. 
Notice that the conventional cone liner model coupling has roughly twice 
the pressure drop of the unskirted bellows while the new configuration had 
about one-half the pressure loss of the bare bellows. 

Another type of liner which appears to have merit is a screen or 
wire cloth liner. Because of its flexibility, a cylinder of wire cloth can 
be fitted directly against the inside of the bellows, allowing some angulation. 
One test has been performed with this type of liner and the results are 
shown in Figure 8. Note that a small increase in pressure loss was realized 
over a bare bellows. It is anticipated that use of a finer mesh wire cloth 
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Figure 61. Sketch Of Bellows Convolutions Giving Unfavorable And 

Favorable Pressure Loss Characteristics 
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Figure 62. Lw Pressure Drop Coupling Design #1 





than used in this first lest will result in an actual pressure loss reduction. 
Possibly of most value, however, is the fact that this type of liner suppresses 
How- induced vibrations; this fact has been verified by testing. Although we 
experienced no problems in the laboratory, the reader will likely realize 
that there are potential problems with a wire cloth liner, should it break 
loose and be carried downstream in the duct system. Attempts to use this 
concept should be thoroughly verified with testing. 

Two other types of liners were tried, also; a longitudinal wire 
screen leaf type, and a spiral wound type using either wire screen or 
metal shim for the liner material. Figures 64 and 65 show examples of 
the two types. In general, all liners suppressed flow induced bellows 
vibrations. For example, peak flow induced dynamic stress over the test 
flow range, in coupling PN 08046 was less than 150 psi with a longitudinal 
leaf screen liner compared with about 10, 000 psi with no liner. 

The pressure loss characteristics of the various liner types were 
quite different. The spiral wound metal shim liner (Figure 66) exhibited 
the best pressure loss characteristics but still needs some development 
work. The problem with the spiral wound metal shim liner was that for 
large bend angles the winding tended to be displaced too far axially, causing 
separation between adjacent windings. If this problem can be cured, this 
type of liner has good possibilities, The spiral wound screen liner (test 
fa not shown) had higher pressure loss than the metal shim liner. 
i igure 67 shows test results for the longitudinal leaf screen liner. In 
general, the pressure loss characteristics at high flow rate are better than 
the unlined bellows for the straight bellows case. For large bend angles, 
however, the longitudinal leaf screen liner had slightly higher pressure 
loss than the unlined bellows at the same angle, 

IV, 4 how Pressure Loss Duct Bend 

Introdu lion 


Duct elbows are known to liave a large pressure loss relative to 
the normal friction loss in straight duct sections. The reason for these 
large losses can be attributed to secondary flows caused by centrifugal 
forces and centritugal force gradients acting on the fluid particles in the 
bend . 


Two experimental elbows have been constructed and Los ted in 
order to evaluate a new method of reducing pressure losses and improv- 
ing velocity profiles in bends, Doth elbows were built with square entrance 
and exit sections to simplify construction, and were identical in every 
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respect, except that the cross-sectional shapes were different in 
the bend region. One elbow was conventional in design and had a 
uniform square cross section. Figure 68 illustrates the new duct 
design. Note that this new design has an increased area on the inside 
of the bend but a decreased area on the outside of the bend; the reason 
for this shape will be discussed below. It was found from testing that 
the new design had roughly 40 percent leas pressure loss than the 
conventional design. It is anticipated that a greatex’ reduction in loss 
could be achieved with more detailed changes in the cross sectional 
shape . 


Low Loss Bend Concept 

Consider Figure 69, which illustrates the flow, of an ideal (no 
viscosity) fluid through a bend. The predicted velocity profiles 
immediately upstream and downstream of the bend are uniform. 

Typical profiles for a real fluid flowing through a bend are shown in 
Figure 69. Here, because of the centrifugal force gradient which 
exists at the start of the turn, the high energy fluid tends to migrate 
to the outside of the elbow while the low energy fluid goes to the inside 
of the elbow; such a fluid migration constitutes a secondary flow. This 
mixing plus the mixing which occurs downstream as the fluid develops 
a normal velocity profile in the downstream duct causes a pressure loss 
which is nominally about an order of magnitude greater than the normal 
friction or viscous losses in a straight duct. 

If on' could, by some means, effectively reduce the driving force 
for this secondary flow without reducing the mean fluid velocity in the 
channel, then the overall elbow pressure loss could be reduced, To 
accomplish this reduction in secondary flow it would be necessary to 
decelerate the fluid near the suction wall and accelerate the fluid near 
the pressure, or outer, wall. Once this were accomplished, there 
should be two results: First, the overall pressure Differential across 
the bend should be decreased and, consequently, low energy fluid migra- 
tion to the inner wall reduced. Secondly, the centrifugal force gradient 
should be reduced sp that the migration of high energy fluid to the outer 
wall is suppressed. 

The manner in which the duct cross .sectional shape tmould be 
ch&ngod to bring about these local accelerations and decelerations of the 
fluid can be estimated analytically. 







Neglecting viscous effects, the velocity profile at the start of a 
bend will tend to be of the form (see Figure 69a) 

k l 

v(r) = ~ 

' r 


where 

_ ( r 2~ r l ) V 
^ “ 4n(r 2 /r 1 ) 

Here, v(r) is the local fluid velocity as a function of radius r measured 
from the center of curvature of the bend, r 2 is the outer wall radius, 
rj the inner wall radius and V the mean fluid velocity in the duct. 

To suppress the secondary ftows it has already been argued that 
the centrifugal force gradient should be made equal to zero (constant 
centrifugal force across the bend); this means that the corrected duct 
local velocity v'(r) should be of the form 



where k 2 is a constant. Therefore, the bend duct height must be 
corrected so that the local velocity in the bend tends to a value 

v'(r) = (k 2 r) 0 ' 5 


\ 

\ 


\ 


or 


v'{r) = 


v(r) k l 

C(r) = rC(r) 


The quantity C(r) is a duct correction factor and has a value equal to 


C(r) = kik 2 “°‘ 5 



5 


To define k 2 , we will arbitrarily spcify that the corrected duct have 
everywhere a centrifugal force equal to the ideal centerline centrifugal 
force in the original duct; in other words 


k 2 



(r 2 +ri ) 3 
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and, therefore 


( r 2 +ri \ 1. 5 
C{r) = 0. 354 ) 


To correct a given duct, the height across the ’-end must be 
modified by the quantity C(r) , In other words, if h(r) is the 
uncorrected duct height as a function of radius, then the corrected 
duct height h'(r) is 


h'(r) = C(r) h(r) 


This type of correction can be applied equally well to square, rectangular, 
circular, etc., ducts. 

Experimental Results 

Two experimental elbows have been constructed and tested to 
verify the concept described in the foregoing. Both elbow shad identical 
cross-sectional areas, and inner and outer wall radii; the only difference 
was that one elbow had a corrected cross-sectional shape in accordance 
with the above. Figure 68 illustrates the corrected elbow. Both test 
items were installed in a water flow system and overall pressure loss 
was monitored. On both elbows, the pressure loss was measured between 
two points several duct widths upstream and downstream from the bend; 
Figure 70 shows the resulting overall pressure loss for both test items. 

It is readily seen that the new design has a lower overall loss, hence, is 
more efficient. 

Figure 71 shows plots of the cross-channel pressure envelopes for 
the two experimental elbows. Note the reduction in the peak pressure 
differential for the new design. This reduction in cross -channel pressure 

dib utial means that the driving force for the fluid migration, discussed 

earlier, has been reduced, and, therefore, a reduction in overall pressure 
loss is expected. 

Conclusion 

The tests performed on the two elb >ws verified that it is possible to 
reduce pressure losses in duct elbows by properly designing she duct cross 
section in the bend region so as. to minimise secondary flows and, hence, 
flow profile distortion. This concept can be applied to a duct with any 
basic shape. 
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Figure 70. Comparison Of Conventional And Low loss Bend Design # 1 
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V, CONCLUSIONS AND RECOMMENDATIONS 


V. 1 Conclusions 


A numoer of conclusions have been derived from the study to 
date; these are listed below: 

Bellows Flow-Induced Vibrations 


(1) Flow-induced vibration of bellows is caused by a fluid- 
elastic instability; specifically it is caused by fluid vortex 
shedding from the bellows convolutions which couples with 
the convolution vibrations. In the absence of convolute 
motion there is no well defined large scale vortex shedding. 

(2) The frequency of the vortex shedding can be related to the 
fluid velocity and convolute dimensions through the use of 
a Strouhal number given by 

S = f s g/V 

where f s is the shedding frequency, V the velocity, and 
a the convolute tip width. This Strouhal number appears 
to be a weak function of only one parameter which is \/o , 
the ratio of pitch to tip width. 

(3) The fluid pressure force which is exerted on a single 
bellows convolute can be expressed in the form 

F = Cp Ap (1/2 pV 2 ) 

where Cp is a vortex force coefficient, Ap is the height 
area of the convolute, p is fluid density, and V velocity. 

The data obtained for Op from a number of bellows tests 
shows it to be a strong function of the parameter X / a . 

(4) It appears that, because of a vortex reinforcement phenomena, 
there are "optimum" values of \ /a which result in local 
minimum or maximum values of Cp. Therefore, there might 
be some advantage in designing the convolute geometry to 
give a minimum Cp value. In general, "open" convolution 
designs are better than "closed" designs. 
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The yespo.jse of a given bellows to flow excitation can 
be predicted with fair accuracy with a method described 
in this report. Strictly speaking, and because no support- 
ing data is presently available, this method is presently 
limited to bellows flowing liquids and having a size of 
about 2-inches ID. We believe, however, the method can 
safely be extended to include gases and larger bellows. One 
important re&triction is that no radial acoustic resonance 
be present to cause amplification of the forces. Future 
work is intended to extend the validity of this bellows flow 
excitation method. 

(6) A very important factor in bellows vibrations is the value 
of the dynamic amplification factor Q. Test results show 
that Q is dependent on the bellows specific spring rate 
(K^N c /D m ), the internal fluid media, and the number of 
plies. A preliminary set of universal Q valuj curves has 
been prepared in terms of known bellows and flow parameters, 
and can be used for present best-eatimate purposes. 

(7) A special parameter, called the "Stress Indicator" has been 
derived and is felt to be of value in determining the relative 
vibration severity of various bellows. All information 
required to make "Stress Indicator" calculations is contained 
in the report. 

(8) A preliminary curve giving fatigue failure information has 
been compiled and can be used to estimate the possibility of 
failure of a given bellows. 

(9) Various kinds of liners can be very effective in suppressing 
bellows fUw-induced vibrations. Care must be exercised, 
however, to keep from causing an increase in pressure loss 
because of the liner. 


(10) Various kinds of externahdamping devices can be effected 

in reducing bellows vibrations. Care must be taken to ensure 
that damping is added to the bellows . 

(11) Rather large errors can be expected when using existing 
spring rate and stress (in terms of deflection) calculation 
methods, largely because of an inaccurate knowledge of true 
convolute dimensional data. 
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(12) The use of an accelerometer for monitoring bellows 
flow -induced vibrations can give misleading results 

as to the moat severs modes o£ vibration, particularly 
where acoustic resonances occur in the system* 

Bellows and Elbow Pressure Loss 

(13) Only one existing bellows pressure loss correlation 
method (see Figure 53) is considered valid for general 
bellows geometry, and this method can lead to rather 
large errors in certain cases, 

/ 

/ 

(14) Careful attention to convolution geometry design can 
result in a 50 percent or more reduction in pressure 
loss when comparing a good design and a poor design. 
Typical good and bad designs were shown in Figure 61. 

In general, a good design is one with small, open convo- 
lutions. Large, close pitch type convolutions generally 
give large losses. 

(15) A bend with changes in cross section such as described 
in this report can result in a substantial reduction in 
pressure loss compared with a bend of constant geometry 
cross section. 

V. 2 Recommendations 

At this point we have several reconimendations: 

. (1) It is recommended that the bellows flow-induced vibration 
calculation procedures presented in this report be used 
for predicting the adequacy of bellows designs and appli- 
cations. The reader should be aware, however, of the 
stated limitations of the method. In particular, care 
j should be taken in those cases where acoustic resonance 

does, or may, occur. 

(2) It is recommended that the use of external damping devices 
be considered as a practical means of suppressing flow- 
induced vibration levels. 


f 
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(3) 


I 




* 


It is recommanded that the bellows pressure loss 
correlation given in Figure 58 be used for prediction 
purposes. 

(4) It is recommended that more work be done to explore 
flow-induced vibrations of much larger bellows and thus 
validate' the prediction method for this case. Also, a 
careful study should be made o£ the acoustic resonance 
phenomena. 

(5) Extreme care should be exercised in interpreting acceler- 
ometer data which is taken fcr bellows flow-induced 
monitoring purposes. Strain gages should be used where 
feasible. 
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APPENDIX A 

dimensional Data for Test Bellows Cited 
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.24 .095 


2 .17 .24 .140 .006 


3 .15 .24 .140 


1 .22 .17 .090 .013 


Note: These dimensions are for geometric description 

only. Actual dimensions for a particular bellows 
may vary somewhat from these representative 
dimensions. 
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APPENDIX B 


Bellows Spring Rate and Stress-Deflection 


BELLOWS SPRING RATE AND STRESS-DEFLECTION 


Spring rate and stress -deflection relationships are of interest in 
bellows vibration studies as well as in quasi-static situations. However, ' 
agreement between experimental and analytical methods is quite poor 
unless a great deal of care is exercised in obtaining the bellows dimensions 
necessary for calculations. Also, the fact that material thickness is not 
constant in a formed bellows greatly complicates any analytical method, 
and necessitates a computer approach or the use of empirical correction 
factors for accurate results. 

Some fairly simple analytical methods are available which can be. 
quite useful if accuracy requirements are not too stringent. 

Axial Spring Rate 

Probably the most successful analytical method for determining 
spring rate is the one developed by Battelle Memorial Institute. The' 
computer program NONLIN which was used in this method provides for 
nonlinear elastic deformations as well as for linear deformation. The 
program is completely described in Battelle's Technical Report No. 
AFRPL-TR -68-22, 

It was found by Battelle that the dimensional data entered into the 
program is of critical importance if reasonable accuracy is expected. In 
fact, the destructive inspection method used in obtaining bellows dimensions 
is considered absolutely necessary. The bellows were encapsulated in a 
potting compound, sectioned along the longitudinal axis, polished, and then 
measured with a machinists microscope. This was done to one sample from 
each manufacturing lot. The remaining bellows in the lot were then-used 
to obtain the experimental results. 

For formed bellows, this technique resulted in calculated values 
which were consistently within 10% 6£ the experimental values. For 
welded bellows the results were less accurate but still quite good. 

Non-computer techniques are usually similar to or based on equa- 
tions developed by F. Salzmann (10), The Salzmann method is outlined in 
Figures B. 1 and B.2. These equations were developed for bellows with 
constant wall thickness and U-shaped convolutions. The results are good 
when this type of bellows is analyzed. This method should be especially 
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Figure B. 1. Bellows Axial Spring Hate 
(From Reference 17} 
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useful with deposited bellows or accurately formed bellows where 
material thinning is not excessive. The importance of accurate 
dimensional data also applies to this method. 

Other similar equations for unusual convolution shapes are 
given in an article “Bellows Spring Rate, “ by J, D. Matheny in Machine 
Design , January 1962. 

The Salzmann method given above is taken from the SAE public 
cation ARP-735. 

When it is impractical to obtain accurate bellows measurements, 
and when approximate spring rate is sufficient, the simplified equation 
shown below is very convenient: 

Ka = D E — P— A\ 3 
A m N c \F/ 

The nomenclature here is the same as in Figure B. 1. 

Table B. I gives a general idea of what can be expected from the 
Salzmann and the simplified methods. 


Convolution Shape 

S 

U 

U 

U 

S 

N c 

9 

8 

7 

8 

21 

N p 

1 

1 

1 

1 

1 

D? 

2.00 

1.5 

1.5 

1.5 

1.5 

D m 

2. 16 

1.74 

1. 86 

1.75 

1.71 

t 

0.006 

0.013 

0. 013 

0.013 

0. 009 

h 

0. 162 

0.256 

0. 345 

0.267 ! 

0.201 

Rate (Test) 

282 

1000 

544 

877 

251 

Rate (Salzmann) 

378 

980 

65 0 

870 

300 

Rate (Simplified) 

364 

855 

42 8 

755 

224 


TABLE B. 1 Comparison of Calculated and Experimental 

Spring Rate 

If very accurate results are required, the spring rate should be determined 
experimentally for each individual bellows. 
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Str e s s -D eflection 


In general, the information given in the discussion of spring 
rate also applies to stress -deflection relationships. Battelle Memorial 
Institute has obtained very good agreement between theoretical and 
experimental methods through the use of a computer program. This is 
reported in the Battelle Technical Report No. AFRPL.-TR-&822, The 
Salzmann method for stress calculation is given in References 10 and 17. 

For the bellows vibration studies in this report, strain gages 
were used mainly as monitors for vibration frequency, displacement, 
and relative stress level. However, a few tests were performed which 
compared experimentally determined stress with that given by the 
Salzmann equation. The experimental values obtained were in variance 
by from about 10 to 30 percent compared with the calculated values. 

In addition to the work of Salzmann and Battelle, other literature 
of potential, value is that of Clark (18), Laupa and Weil (19) and Turner 
and Ford (20). 
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APPENDIX C 


Prediction of Bellows Added Fluid Mass 
for Out-of-Phase Mode 



PREDICTION OF BELLOWS ADDED FLUID MASS 
FOR OUT-OF-PHASE MODE 


Consider the idealized situation shown in Figure C. 1; here we 
have two adjacent convolutions vibrating out of phase. It is assumed 
that the convolution sides remain straight at all times. This motion of 
two convolutions toward one another forces some fluid to move in the 
y direction with a velocity as a function of y given by 


v(y) - 


xy 2 

6h 


(C.l) 


From the momentum equation, this fluid motion will cause a pressure 
gradient along the convolution given by 


or 


9£ 

9y 



„2 

dp _P£ 
dy = 6h 


M 


X 


(C. 2) 


(C.3) 


Integration gives (see Figure C.l) 


p(y) = 


PfX 

"36' 


h 2 -^ 


(C.4) 


This is, then, the pressure acting along one side of the convolution 
because of out-of-phase acceleration of two adjacent convolutions. The 
problem now is to compute an effective force which acts on the convolution 
as it vibrates. If we knew precisely the mode shape of vibration of the 
convolution (it has been assumed linear) then this effective force could be 
calculated exactly (as could the fluid velocity distribution); we don't know 
the mode shape, however, so that the effective force must be estimated. 

It is reasonable to assume that it falls somewhere between two limiting 
values. The first is computed by calculating the total force acting over 
the convolution (both sides) which gives 


F 


^ ^ m.Pfk 


i* 

X 


26 


(C . 5) 



on or Fluid Added Mass trreci 
ir Out -of -Phase Mode 

a 

.1 

. « 

iJ 

a 

n 

i’ 

1 

y 

e 

n 

C2 ■' s 

0 

n 



> 


This value of force is believed to be too high. A lower limiting value 
is felt to be obtained by weighing p(y) with y and integrating over 
the convolution, hence 

p(yiydy 

Tr D m Pfh 3 „ (C.6) 

= g" 1 ■ x 

5o 


A more reasonable value is believed to be 

3 


F = 


TrP m Pfh- 
36 


(C.7) 


which gives an effective added mass for the out-of-phase mode of 


m 


o 

7TD m p£h 
— — 


(C.8) 


For two different cases, use of this value of added mass has given good 
results when experimental resonant frequency data was compared with 
calculated data; one was the ring convolution test model discussed 
previously, the other was a flexible hose for which the in -phase and 
out-of-phase modes were both calculated and observed. 


C3 



APPENDIX D 

Excitation of Longitudinal Acoustic Resonance 
in Plow System 



EXCITATION OF LONGITUDINAL ACOUSTIC RESONANCE 

IN FLOW SYSTEM 


When a bellows in a fluid line is experiencing vibrations induced 
by the fluid flow, it is generally noted that strong acoustic longitudinal 
standing waves are excited in the flow system piping. This phenomenon 
is caused by the bellows vibration exciting some fluid system acoustic 
mode(s), or more generally stated, resulting in a coupled bellows -liquid 
mode of the type discussed in the previous section. This coupling effect 
may be important for two reasons; first, it is feasible that the acoustic 
properties of the entire flow system may have a significant effect on the 
bellows flow-induced vibration characteristics and, second, strong 
acoustic resonance in the flow, with the bellows as a source, may 
adversely affect the performance of other components in the system, 
for example, the tur'bopump. 

For the present discussion, the problem will be simplified by 
considering the bellows vibration as uncoupled from the flow system 
acoustics and simply representing a sound source in the system. 

Consider Figure D. 1 which shows an idealized case in which we 
have a ring vibrating in a pipe; this simulates the vibration of a single - 
bellows convolution. From an acoustic standpoint, this vibrating ring 
is a dipole sound source. Without going into the mathematical details 
of the near field sound characteristics of this source, we will assume 
that, for all practical purposes, this vibrating ring produces a fluctuating 
difference in pressure between a point immediately upstream and another 
point immediately downstream of the ring. 

Now, consider Figure D.2, which shows a simplified schematic 
of our flow system with this idealized dipole source located as shown. 
Acoustically, we are assuming the system to consist of a line or duct with 
a dipole source at some intermediate point and with a constant pressure 
reservoir at each end. It is important to note that the pressures p and 
velocities v shown in Figure D.2 denote acoustic or nonsteady disturbance 
quantities, not total quantities. This explains why the acoustic pressure 
is zero in each reservoir, even though the steady-state pressure is not 
zero. The problem now is to calculate the response of the fluid system 
to this dipole source. 

With little loss of accuracy, the fluid lines on either side of the 
source (Lj and L 2 ) will be assumed frictionless and one -dimensional from 
an acoustic point of view. The acoustic impedance of each line (L^ and L^) 
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looking away from the vibrating ring sound source is 


p 2 (iw) 

v^ET) = - ipc tan( " Ij i /c) 


and 


p 3 (iw ) 

—7— = ipc tan(wL 2 /c) 

y 3 (xw) * 


(D.l) 


(D. 2 ) 


where p is fluid density and c speed of sound. These, two relations 
can be combined to give 

p 2 (ioo) -p 2 (iw) = -ipc v 3 (iw) |tan(wLj /c) + tan(wL 2 /c)j (D. 3 ) 


From Equation (D. 3 ) the magnitude and phase of v 3 can be determined for 
a given p2 - P3 = A e iwt (value of sound source). Knowing p 2 (iw), 

P 3 U w ) and v 3 (iw ), then the acoustic pressure and velocity at any position 
X 1 or x 2 along the ducts can be determined from the equations 

p(xj, iw) =p 2 (iw) cos[w(L 1 -x 1 )/c] + ipc v 3 (iw) sir>[w(L 1 -Xj )/c] (D. 4 ) 

v(xj, iw) =v 3 (iw) cos[w(Li-xj)/c]+ p 2 (iw) sin[w(Li~x 1 )/c] (D.S) 

r 

and 

p(x 2 , iw) = p 3 (iw) cos (wx2/c) - ipc v 3 (iw) sin (wx2/c) (D.6) 

v(x2>iw) = v 3 (iw) cos (wx2/c) - p 3 (iw) sin (wx2/c) • (D. 7 ) 


The above relations describe standing waves in the flow system ducts. The 
properties of these standing waves, that is the location of pressure nodes, 
antinodes and possible resonances, depend essentially upon the values of 
Lj and L 2 « 
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Figure D.3 System Acoustic Standing l 
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For example, assume Lj = L 2 = L, and also put X 3 = - X^; 

eliminating P 2 , P 3 and V 3 from Equations (D.4) through (D.7) gives 

(D. 8 ) 
(D.9) 
(D. 10) 

and 

v(x 2 ,iio) =v(x 1 ,iw). (D. 11) 


A e iwt ( 1 

p(xj . iu) = — ~2 ~~ ) cos ( wx 3‘ /c ) -sin(o>X 3 /c) • cot(wL/c)j 

iA e itot t f 

v(xj,iw)= 2 pc "jcosfwx^/c) cot(wL/c) + sin{wx 3 /c)j 


p(x ? , iw)= -p(x 1 ,iw) 


Notice that the system will exhibit a resonance condition when 

2L 

toL/c = ntr, (n is any integer); antiresonance occurs when wL/c = 2 = nTT. 
Figure D.3 shews the system standing wave patterns predicted by Equations 
(D. 8 ) through (D. 11) with uL/c = 5ir/4 as an illustration of a typical situation. 
Note the jump discontinuity in pressure at x = 0 with the value of the jump 
being equal to the strength of the dipole source. Further, note that the 
peak system acoustic pressure does not occur at the source, but rather 
occurs at some point along the line(s); this fact is generally true whenever 
coL/c > tt/2. 

As mentioned previously, the above formulation considered the 
vibrating ring (bellows idealization) to be simply a sound source for the 
flow system, hence, uncoupled from the system acoustic properties. In 
the more general case, the influence of the system upon the bellows 
vibration should be considered. 
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